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TECHNICAL NOTE 3970 

THIN AIRFOIL THEORY BASED ON APPROXIMATE SOLUTION 
OF THE TRANSONIC FLOW EQUATION 
By John R. Spreiter and Alberta Y. Alkane 


SUMMARY 


The present paper describes a method for the approximate solution of 
the nonlinear* equations of transonic small disturbance theory. Although 
the solutions are nonlinear, the analysis is sufficiently simple that 
results are obtained in closed analytic form for a large and significant 
class of nonlifting airfoils. Application to two-dimensional flows with 
free-stream Mach number near 1 leads, for instance, to general expressions 
for the determination of the pressure distribution on an airfoil of spec- 
ified geometry and for the shape of an airfoil having a prescribed pressure 
distribution and gives, furthermore, the correct variation of pressure 
with Mach number at Mach number 1. For flows that are subsonic everywhere, 
the method yields a pressure-correction formula that is more accurate than 
the Prandtl-Glauert rule and compares favorably with existing higher 
approximations. For flows that are supersonic everywhere, the method 
yields the equivalent, in transonic approximation, of simple wave theory. 
Results obtained by application of these general expressions are shown to 
correspond closely to existing solutions and to experimental data for a 
wide variety of airfoils. 


INTRODUCTION 


The difficulty of solving the nonlinear equations of motion of a 
compressible inviscid gas has led to widespread use of approximate methods 
in the practical solution of the problems of airfoil theory. The simplest 
and most versatile approximate method is that based on a complete lineari- 
zation of the equations and stems from the pioneering work of Muhk, Prandtl, 
Glauert, Ackeret, and others (see refs. 1 and 2 for a resume). Although 
this linear theory of compressible flow has been extensively developed in 
recent years and is widely used in aeronautical applications, it has two 
limitations that are of significance in the present discussion. First, 
linearized theory gives only a first approximation that is correct for 
airfoils of small thickness ratio. This limitation is, in some respect, 
of continually diminishing significance as the aeronautical engineer is 
forced to use thin wings and slender bodies to avoid heavy penalties in 
wave drag. If the airfoil is not sufficiently thin, however, corrections 
are necessary and higher order theories have been developed to fill the 
need (see ref. 3 for a resume) . Second, and more important for the 
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present discussion, linearized theory requires, in general, that the Mach 
number he sufficiently removed from unity that the flow is either purely 
subsonic or purely supersonic. If "both subsonic and supersonic velocities 
occur in different parts of a single flow field, the flow is said to he 
transonic and the results of neither linear theory nor the existing 
higher order theories are, even qualitatively, in agreement with the 
experimentally observed flows. 

Transonic flows have been studied successfully by consideration of a 
simplified, although still nonlinear, theory that was originally conceived 
in an effort to provide a useful first approximation for the pressures and 
forces on thin wings and slender bodies in inviscid flows with free-stream 
Mach number very near unity (see ref. 4 for a short resume) . More recent 
developments described in references 5> 6, r J } 8, and elsewhere have shown 
that the useful range of this theory can be extended to include subsonic 
and supersonic flows if slightly different approximations are employed in 
the derivation of the fundamental equations. Although the resulting theory 
is commonly designated as transonic small disturbance theory, or more 
briefly as transonic flow theory, it is actually a unified theory for sub- 
sonic, supersonic, or transonic flow around thin wings and slender bodies, 
a.nrl is moreover, the silkiest theory proposed to date that is capable of 
yielding reliable results throughout that Mach number range. 

This formulation of transonic flow theory provides a set of equations 
that differs from that of linear theory by the addition of one nonlinear 
term in the differential equation for the perturbation potential and in 
the shock relation. If the flow is purely subsonic or purely supersonic, 
solutions of the equations of transonic flow theory can be sought by appli- 
cation of existing methods for approximating the solutions of the exact 
equations of compressible inviscid flow. If the flow Is transonic, how- 
ever, the results obtained by application of these methods are at wide 
variance with those observed experimentally and it is necessary to devise 
new and appropriate methods of solution. Although methods of the succes- 
sive approximation type have recently been developed that can be applied 
to transonic flows (e.g., refs. 9 end 10), the principal method that has 
been employed in the theoretical analysis of such flows involves the use 
of the hodograph transformation by means of which the nonlinear equation 
for the perturbation potential is transformed into a linear differential 
equation of mixed elliptic-hyperbolic type, the Tricomi equation. Although 
the resulting boundary-value problem is still very difficult to solve, 
this method has been applied with considerable success in the study of 
transonic flow around wedge and flat-plate airfoils and a number of spec- 
ific results have been given in recent years by Guderley and Yoshlhara, 
Vincent! and Wagoner, Cole, and others (see ref. 11 for a resume). Exten- 
sion of this method to permit calculation of transonic flows around arbi- 
trary airfoils with curved boundaries appears, however, to be a difficult 
task. 


The present analysis is based on a novel method of approximation that 
avoids most of the difficulties of existing procedures while still preserv- 
ing much of the nonlinear effects in the solution. Sufficient simplicity 
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is gained "by restricting attention to surface pressures and to flows that 
are either purely subsonic, purely supersonic, or have a free- stream Mach 
number near 1, that results can be obtained in closed analytical form for 
both the direct problem of calculating the pressure distribution on an air- 
foil of given shape, and the inverse problem of calculating the shape of 
an airfoil associated with a given pressure distribution. Inasmuch as the 
magnitude of the errors introduced by use of the approx im ation procedures 
is not evaluated in all cases by mathematical considerations, the useful- 
ness and accuracy of the results are demonstrated by the calculation of 
the pressure distribution and drag for many different airfoil shapes and 
by comparison with existing theoretical and experimental results. 

Of the theoretical results available for comparison, only two are 
exact. They are the simple-wave solution for supersonic flows without 
shock waves, and the variation of pressure with Mach number at Mach num- 
ber 1. The present method yields both of these results exactly within the 
framework of transonic small disturbance theory. 

Although the existing results mentioned above for wedge airfoils at 
Mach number 1 contain certain approximations beyond those implicit in the 
use of the equations of transonic flow theory, the influence of these 
approximations appears to be minor and the results are generally consid- 
ered to be very nearly exact solutions of these equations. The present 
method produces results for this case that are in substantial agreement 
with these previous theoretical results. In contrast to the holograph 
methods, the necessary steps are sufficiently simple, moreover, that 
results can also be obtained for sonic and near sonic flow around arbitrary 
airfoils with curved boundaries. Since previous theoretical information 
for such cases is meager, comparisons are made with a large number of 
experimental results. In general, the theoretical results found by appli- 
cation of the present method lie within the range of experimental scatter 
of the data. 

In the subsonic range, no exact solutions are available for flow 
around a thin airfoil. Comparisons are made, therefore, with pressure 
correction formulas, such as that of Karman-Tsien, and with higher 
approximations obtained by iteration methods. 

A simple heuristic account of the general method and extensive dis- 
cussion of the results are contained in the main text. Additional details 
concerning the underlying basis for the general procedures are contained 
in the Appendix. 


PRINCIPAL SYMBOLS 


a 


speed of sound 

speed of sound in the free stream 
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U E 

U H 

Up 
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Uoo 

local Mach number 
free- stream Mach number 

exponent in the relations for airfoil orainates given by equa- 
tions (56) ana (6l) 

static pressure 

free- stream static pressure 

resultant local velocity 

maximum thickness of profile 

free-stream velocity 

perturbation velocity components parallel to x ana z axes, 
respectively 

value of u obtal nea by solution of equation (21) 

value of u obtainea by solution of equation (8) 

value of u obtainea by solution of equation (39) 

value of x at which parabolic ana hyperbolic solutions are 

joinea 
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x, z Cartesian coordinates where x extends in the direction of the 

free-stream velocity 

x* value of x at which the local velocity is sonic 
Z ordinates of the upper surface of the airfoil 

T gamma function 

7 ratio of specific heats , for air 7 * 1.4 

9 semiapex angle of wedge airfoil 


A 

-r\y ^ r pm "pp cry* or] ptvP. _ 

ao p 



d(x/ c ) 

?*E 

1 - l^ 2 - ku 



M 2 - 1 + ku 
00 


Xp 

k^ 

Sx 


? co 

M 2 - 1 
00 

[1^ 2 (7+1)t] 2 / 3 


p 

00 

free-stream. density of air 

T 

thickness ratio, t/c 



9 perturbation velocity potential 

Subscripts 


cr values associated with critical Mach number ■ 

i values associated with incompressible flow or with 

L values given by linearized compressible flow theory 

M *1 values associated with M =1 

CO 00 

Zjaax values associated with maximum ordinate of airfoil 
| =0 values associated with | = 0 , or with M =1 

00 3 CO 7 CO 
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FUHDAMEHTAL EQUATIONS AND BOUNDARY CONDITIONS 


Consider the steady flow of an Inviscid compressible gas past an 
arbitrary thin symmetrical nonlifting airfoil, and introduce Cartesian 
coordinates x and z with the x axis parallel to the direction of the 

free- stream, as illustrated in sketch (a) . 

Let the free-stream velocity and density 
be Uqo and p M , the perturbation potential 
be cp, and the perturbation velocity compo- 
nents parallel to the x and z axes be 
or u, and cp z , or w, where the subscript"' 
indicates differentiation. The boundary con- 
ditions require that the perturbation veloc- 
ities vanish at infinity, and that the flow 
be tangential to the wing surface. The first 
is constant at infinity. The latter condition 



V 


Sketch (a) 


condition indicates that cp 
can be approximated for thin wings by 


v 


(<P ) 

z=o 


U„ 


dZ 

dx 


( 1 ) 


where Z represents the ordinates of the airfoil upper surface . The 
pressure coefficient Cp is likewise approximated to first order by 




P-P« 





( 2 ) 


These relations are familiar from linear theory, but apply equally for 
transonic thin airfoil theory. The differential equation for cp is not 
the same as in linear theory, however, but is 


<P»-° ^ 

where is the Mach number of the undisturbed flow and 7 is the ratio 

of specific heats (l.l for air). It is useful to note that the coeffi- 
cient of cPxx corresponds, in the present approximation, to 1-M 2 
where M is the local Mach number. 


Knowledge of methods for obtaining solutions of equation (3) is 
meager, not only because the equation is nonlinear, but because It can 
change type (elliptic, hyperbolic) , depending on the value of 1^, and <P X . 
This change of type is an essential feature of transonic flow, since 
subsonic flows are represented by elliptic equations and supersonic flows 
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by hyperbolic equations. If both types of flow occur in a single flow 
field, it is apparent that the differential equation must change type. 
In the present ease, the type of the equation is recognized by the sign 
of the total coefficient of cp^, as follows: 


He 2 - H , 2 


> 0 elliptic (subsonic) 

< 0 hyperbolic (supersonic) 


(4) 


An important quantity in the discussion of compressible flows is the 
critical pressure coefficient C-n associated with the local occurrence 

■t'er 

of sonic velocity. The appropriate relation is found by combination of 
equation (2) and the relation obtained by equating the coefficient of cp 
in equation (3) to zero, and is ** 



2(1 - O 

H» 2 (7 + 1) 


(5) 


In transonic and supersonic flows, it is also necessary, in general, 
to provide an additional equation for the discontinuous changes in veloc- 
ity that occur at shock surfaces. The necessary equations, when simpli- 
fied to the form consistent with the approximations of transonic flow 
theory, reduce to 


i-^ 2 -^c 2 ( 2S±3l ) (u a -u t ) 2 + (w^) 2 - 0 


( 6 ) 




where the subscripts a and b refer to the values on the two sides of 
the shock surface. With the exception of the Appendix and minor refer- 
ences in the main text, equation (6) is not employed explicitly in the 
following analysis and discussion because attention is confined to 
(a) purely subsonic flows in which no shocks occur, (b) purely supersonic 
flows in which shock waves can be approximated with good accuracy by isen- 
troplc compressions, and (c) flows with free-stream Mach number near 1 in 
which the shock waves are situated either downstream or far upstream of 
the airfoil. 


The remainder of the present paper is concerned with the approximate 
solution of the precedi n g equations and with comparison of the results 
obtained in specific applications with existing theoretical and experi- 
mental results. Purely supersonic flows are discussed first because the 
method of approximation yields the exact equivalent, in transonic approx- 
imation, of simple wave theory. Purely subsonic flows are discussed next 
because of the close relationship between the results for this and the 
preceding case. Flows with free-stream Mach number near 1 are treated 
last. 
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SUPERSONIC FLOWS 

Approximate Solution of Equations 


It is convenient in the analysis of supersonic flows to introduce 
the symbol Ag as an abbreviation for the negative of the coefficient 
of 




M 


1 + M oo 2 -TT <P X * ^o 2 - 1 + ^ >° 


(7) 


and rewrite equation ( 3 ) in the form: 


Tfl^xx + ?zz = 0 


( 8 ) 


It is now assumed that A H is neither zero nor infinite and that it varies 
sufficiently slowly that its derivatives can be disregarded so that it can 
be considered, temporarily, as a constant. At this stage, the problem is 
equivalent to that encountered in linearized supersonic airfoil theory 
(it is identical if Ag is replaced by M^-l) and the solution Ug at 
the airfoil surface is 


u H 


U oc dz 


Differentiation yields 


du H 

dx 


Uqq d a Z 


(9) 


( 10 ) 


If, now, W^-l + ku is restored in place of Ag so that, in effect, 
the local value for Ah is used at each point and the subscript H 
on u is dropped, equation (10) becomes 


Jkg. gfg (11) 

^ dx2 

It is immediately apparent that a certain degree of arbitrariness is dis- 
played in the preceding steps and that different results will be obtained 
depending, for instance, or whether M^-l + ku is substituted for Ah 
in equation (10) as above, or in equation ( 9 ), or in other equations 
obtained by further differentiation or integration of equation ( 9 ) • It 
is shown in the Appendix, however, that the error involved in the preced- 
ing steps can be assessed exactly by examination of the remainder terms 
that have been omitted in writing equations ( 9 ) and (10). The advisabil- 
ity of using equation (ll) is assured by the fact that the error is .shown 
to vanish, in the absence of shock waves, if Ah is replaced by l^-l+ku 
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in equation (10) , tut not in equation (9) • This conclusion "becomes imme- 
diately evident, furthermore, upon recognition of the fact that equa- 
tion (ll) is the counterpart, in transonic small disturbance theory, of 
a fundamental differential equation that occurs in the analysis of Prandtl- 
Meyer and simple wave flows (see, e.g., ref. 12, p. 87 or ref. 13, pp. 190 
and 212 ) . Equation (ll) is a nonlinear ordinary differential equation 
for u that can be solved easily be separation of variables. The result 
is 


4 


| - (V-l+lni ) 372 =. -U M + C (12) 

3 h dx 

where C is a constant of integration. In applications of equation (12) 
to flows that are supersonic everywhere, perhaps the most logical method 
for the evaluation of this constant is to use the expression between u 
and dZ/dx provided at the leading edge by the transonic approximation 
to the shock relation, that is by equation ( 6 ) with and w a equated 
to zero, uj, to (u) x=0 , and w^, to U co (dZ/dx) x=0 . The result given by 
equation (12) with C evaluated in this way corresponds, to the degree 
of approximation afforded by -use of transonic small disturbance theory, 
to shock-expansion theory. An alternative procedure that leads to a 
somewhat simpler result possessing very nearly equal accuracy is to eval- 
uate C by use of the result indicated by equation ( 9 ) that u = 0 
where dZ/dx * 0 for any nonsingular Ag; thus 


C 


_ 2 _ 

3h 


(Hjo 2 - 1 ) 


3/2 


(13) 


whence 


u 


1 

k 


r ,2/3 

+ (W-l) 37 * - I 


(14) 


The corresponding relation for the pressure coefficient Cp is obtained 
by combination of equations (2) and (l4) , and Is 


°P - 


r 


VK7+1) 


(H» 2 -i) “ (i^ 2 -i ) 3/2 - \ n» 2 (r+ 1) f§ 


,2/3 


L 


dx 


(15) 


It should be noted that the restriction to supersonic flow imposed in the 
evaluation of C and in the inequality of equation ( 7 ) requires that 
equation ( 15 ) is to be applied only to cases for which the quantity in 
the square brackets, that Is, [(Moo 2 -!) 3 / 2 - ( 3 / 2 )M O 0 2 ( 7 +l) (dZ/dx) ] , is 
positive. 


* 
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Comparison With Existing Higher Approximations 


Equation (15) is recognized, by comparison with equation (3-15) of 
reference 3* page 387* 1 as the precise equivalent, in the transonic small 

_ c disturbance approximation. 


— Mco-I' 

Mco-U-S, 




Eq. (15) 

Exact simple — 

I I wave theory 


Sketch (b) 


of simple wave theory for 
the surface pressure on an 
airfoil in supersonic flow. 
Exact simple wave theory is 
known, moreover, to be per- 
fectly adequate for all 
practical purposes up to a 
Mach number of 3, which is 
considerably in excess of 
the present range of inter- 
est. Within this Mach num- 
ber range, the results 
obtained by use of simple 
wave theory are almost iden- 
tical with those obtained by 
use of shock- expans ion 
theory. Comparisons of the 
variations of Cp with dZ/dx 
indicated by exact simple 
wave theory and by equa- 
tion (15) are shown in 
sketch (b) for several Mach 
numbers from 1 to 2. As 
might be anticipated, the 
two sets of results are in 
close agreement for Mach 
numbers near 1, and differ 
by an increasing amount 
with increasing Mach number. 


Although the necessary calculations are very easy to accomplish in 
any given case, simple wave theory is not always used in actual practice. 
Many calculations are based on linear theory or Busemann's second-order 


■’•Comparison discloses that the quantity ML 2 (7+1) that appears in equa- 
tion (15) is represented by 7+1 in equation (3-15) of reference 3- The 
difference is associated with a corresponding difference in the coeffi- 
cient k of the nonlinear term of equation (3) . Although the two coeffi- 
cients are identical at ^=1, and might appear to be equally consistent 
with the other assumptions of transonic flow theory, it has been shown in 
references 5, 6, 7, 8, and elsewhere that the approximation obtained by use 
°f M co 2 (7 + i) is much the better of the two for Mach numbers other than 1. 
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Sketch, (c) 

theory. Consequently, an additional set of graphs is shown in sketch (c) 
in which the curves of sketch (h) are repeated together with the corre- 
sponding curves calculated h y use of first- and second-order theory. No 
comparisons are shown for 1^=1 because the latter theories Indi cate 
infinite pressures. It can be seen that equation ( 15 ) furnishes a better 
appro ximat ion than linear theory throughout the entire range of variables 
shown on sketch (c) and a better approximation than second-order theory 
for Mach numbers less than about l.^t-. It can be seen that second-order 
theory furnishes a very poor approximation for Cp at Mach numbers 
approaching unity. 


In order to explore this behavior further, two additional curves 
labeled "third order” and "fourth order,” calculated using the formulas 
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of references ih and 15 j 2 are included on the graph, of sketch (c), even 
though they must he interpreted in a somewhat more restricted sense than 
the other curves. To he more precise, the third-order curve is restricted 
to airfoils for which dZ/dx is zero at the leading edge, and the fourth- 
order curve to airfoils for which both dZ/dx and d £ Z/dx 2 are zero there. 
It is clear from this sketch that the accuracy of second-order theory at 
Mach numbers near unity is not improved by addition of higher order terms. 
The explanation resides in the fact that the larger values of |dZ/dx| 
shown on the graphs of sketch (c) exceed the radius of convergence of the 
power series expansion for Cp for all but the highest Mach number shown. 
With the noted restrictions on the leading edge, the higher order results 
of sketch (c) are equivalent to the first few terms of a power series 
expansion, in terms of dZ/dx, of the expression for Cp indicated by 
exact simple wave theory. The radius of convergence of the series depends, 
of course, on the Mach number and is given by the value of |dZ/dx) asso- 
ciated with the occurrence of sonic flow or, in terms of the curves shown 
on sketches (b) and (c), with the termination of the left end of the exact 
curve. The failure of higher order theories at negative dZ/dx is thus 
of purely mathematical origin and has no direct physical significance. 


Additional Properties of Approximate Solution 


Equation (15) has some additional interesting properties worth noting. 
Of the two major components of the right-hand side, the first is recognized 

upon comparison with equation (5) as the expression for CL. . Since the 

■^cr 

remaining term is zero when CL) = CL, > it follows that the expression for 

-cr 

the critical value for dZ/dx associated with the occurrence of sonic 
velocity at a given Mach number M^ is 


/dZ\ 2(^-1 ) 3/g 

W/ cr * 3Hx, 2 (7+1) 


(16) 


It follows, furthermore, that a curve representing the variation of Cp 
with ^ for a given dZ/dx, and hence a given point on the airfoil, 
approaches infinite slope as Cp approaches Cp . 


An alternative form for equation (15) that is useful for some pur- 
poses is the following which expresses Cp in terms of the linear-theory 
solution Crv. rather than dZ/dx. 


2(H/-1) , L 3 M» 2 (7+l) „ 1 2/3 ] 

1^ 2 (7+1) L £ M* 2 -! ^ J J 


(17) 


Attention of those who refer to reference 15 is called to the fact 
that the first term appearing in the fourth-order coefficient a* of 
equation (29) should be 2/3 rather than 1/3. This term is written 
correctly in the numerical example given in equation (135) » 
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where 


n , 2 iZ 

■^L l v 2 _-, dx 


V^ 2 - 1 

This relation can "be written in somewhat more concise form if expressed 
in terms of the transonic similarity parameters Cp and 1^, thus 

, _ \ a / 3 ~| 

l-'l'T” ) (18) 

4 S»/ J 


- 2 i. 


where 


_ _ [lC(7+^)] 1/a 

°P “ Zz ? 3 °P 


^o g -l 

6 “ " [H/Cr+DU* 73 


nnfi t refers to the thickness ratio. Critical values for Cp and Cp^ 
corresponding to the local occurrence of sonic velocity are easily recog- 
nized to "be the following: 


°p =r - 2! « > 


a - 4 t 

°Pl " Z 5 °° 

Jj cr J 


(19) 


SUBSONIC FLOWS 

Approximate Solution of Equations 


The procedure described in the preceding section will now be applied 
to the analysis of subsonic flows. Thus, introduce the symbol Ag as an 
abbreviation for the coefficient of cp^ 


A e * l-^-Moo 2 ££ 9 X = l-M O0 2 -ku>0 
w oo 


( 20 ) 


and rewrite equation (3) as follows: 

VPxx + 'Pzz - 0 (21) 

If it is again assumed that Ag is neither zero nor infinite and that it 
varies sufficiently slowly that its derivatives can be disregarded, the 
problem is equiv alent to that encountered in linearized subsonic airfoil 
theory and the solution ug at the airfoil surface is 


^5 


U a 


7t -fkft 


r 


dZ/d| 


d| 


ui 



x-g 


( 22 ) 
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where the subscript i refers to the values for M^-O. Differentiation 
yields 


du B 1 du i 

,/*E ^ 


( 23 ) 


If, in the same manner as described for the supersonic case, l-M^-ku 
Is restored in place of Ag so that, in effect, the local value for A® 
is used at each point, and the subscript E on u is dropped, equation (23) 
becomes 


du = 1 dui 


( 24 ) 


As in the previous discussion of supersonic flows, the error terms are 
omitted in writing the preceding relations, but are included in a more 
complete presentation of the equations given in the Appendix. Once again, 
the resulting relation is a nonlinear ordinary differential equation that 
can be solved readily by separation of variables 




ku) 


a/2 


up + C 


( 25 ) 


In applications to flows that are subsonic everywhere, the constant of 
integration C is evaluated by use of the result indicated by equa- 
tion (22) that udO where upoO for any nonsingular Ae> thus 


0 - - 1 d-0 3A 


(26) 


In this way, the following relation is obtained between u and up 

l 2 / 3 


1 

u * — 
k 


l-l^o 2 - 


- | fay. 


( 27 ) 


The corresponding relation for the pressure coefficient Cp is obtained 
by combination of equations (2) and (27) and is 


Cp 


^(7+1) 


(l->4o 2 ) 


2 ' 1 fl K> s f /S + r H» 2 (7+1)C 


- (1-i 


n2/a 


Pi 


( 28 ) 




-2 - - 2 r° dZ/d* 


U n 


X-l 


d| 


where 
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In tiie same way as noted for supersonic flows following equation ( 15 ) , 
the restriction to subsonic flow imposed in the evaluation of C an d by 
the Ine quali ty of equation (20) requires that equation (28) be applied 
only to cases where the quantity in the square bracket is positive. 

This result possesses several simple, but interesting, properties. 
First of all, the leading term of an expansion of equation (28) in a 
series involving ascending powers of Cp^ is precisely the familiar 
Prandtl-Glauert rule of linearized subsonic compressible flow theory 



(29) 


The coefficients of succeeding terms, however, do not agree with those 
given by the method of successive approximation. Next, the first of the 
two major components of the right-hand member of equation (28) is recog- 
nized, Just as in the supersonic case, as the expression for the critical 

pressure coefficient . Since the remaining term is zero when 

-=er 

Cp=Cp , it follows that the expression for the critical value for Cp^ 

associated with the occurrence of sonic velocity at a given Mach num- 
ber 1^, is 

» q-io 3/a 

icr 3 Hx, 2 (7+1) 



It may be noted that this value is Just two -thirds of that obtained by 
use of equation ( 5 ) together with the Prandtl-Glauert rule. It follows, 
furthermore, that a curve illustrating the variation of Cp with f° r 
a given Cp^, and hence a given point on the airfoil surface, approaches 
infinite slope as Cp approaches Cp cr * This latter behavior signifies 

that a power series expansion of the result will only converge for Mach 
numbers less than the critical. Last, the following result is obtained 
if equation (28) is expressed in terms of the subsonic linear theory solu- 
tion Cp^ rather than Cp 


where 


-2(l-M C0 2 ) 

^=o 2 (7+1) 


1- 


" 1+ 1 ^(7+1) 

_ ^ 1-Mco 2 



°PL 


°Pi 

J 1-Moo 2 


(31) 
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Note that the relation between Cp and Cp_ indicated by equation (31) for 
subsonic flows is precisely the same as given by equation (17) for super- 
sonic flows. It follows Immediately 
that the corresponding expression in 
terms of the_transonic similarity 
parameters Cp and | ra given by equa- 
tion (18) applies to subsonic, as well 
as supersonic flows. In order to 
illustrate the nature of the results 
indicated by equation (l8) , a plot of 
the variation of Cp with for 

various Cp^ [1^] is shown in 

sketched) . Although the remarkable 
symmetry about 1,^=0 is a consequence 
of expressing the results in terms of 
the transonic similarity parameters, 
the general symmetry remains, although 
in somewhat distorted form, when Cp 
is plotted as a function of for 

constant Cp^ 1 l-M^j . Such a plot 
is shown in sketch (e) . 


Comparison With Existing 
Higher Approximations 


The remainder of the present sec- 
, > tion on subsonic flows is concerned 

Sketch (>e; with an evaluation of the degree of 

accuracy achieved by use of equation (28) . This discussion is handicapped 
somewhat by the fact that all other theories for subsonic flows around air- 
foils are also approximate and that no exact solutions are known. Perhaps 
the most widely used higher approximation is the Karman-Tsien rule. (See 
refs. 3 and 1 6 for a resume.) Although the traditional derivation of the 
Karman-Tsien rule is based on the hodograph method, it is not without 
interest to observe that the Karman-Tsien rule can be obtained by use of 
the present procedures together with the three assumptions introduced in 
the original derivation. These are: (a) that 7 can be approximated 

by -1 in the expression for the speed of sound, (b) that the perturbation 
velocities are small, and (c) that the Mach number can be considered small 
in the evaluation of additional effects of compressibility beyond those 
indicated by linear theory. The starting point is equation (3) "with the 
coefficient of cp^ replaced by the approximate relation for 1-M 2 , 
where M is the local Mach number, afforded by use of assumptions (a) 
and (b) . The necessary relation can be derived from the energy equation 



as + 


7-1 


7-1 




(32) 
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where 


q. 2 « (Uoo+ u ) 2 + v2 


lay setting 7 * - 1 , whence 


a 2 - a^ 2 - U^ 2 + q. 2 


« 


and 


1-M 2 



l-^ g * l-^o 2 

1-Hx, 2 +H» 2 (q. 2 /Uoo 2 ) l+2M aj 2 (u/U co ) 


(33) 


Note that this approximation does not permit the attainment of M=1 with 
finite q/U w and with I^o different from unity. Substitution of this 
relation for 1-M 2 in place of Ag in equation ( 23 ) and integration 
yields 

^ l+2M O0 2 (u/U co ) - - 7 == ~ + c (3*0 

n/1-Wk, 2 U “ 


The constant of integration is again evaluated by setting u=0 where u-l= 0, 
whence C equals unity. Solution for u and introduction of the rela- 
tion between u and Cp given by equation ( 2 ) yields 


' . -£L- (- CpV ^ 

A-n» a y 4 -/i-h* 2 / Vl-h/ + (i/4)H» s Cp 1 


(35) 


Replacement of M^ 2 in the second term in the denominator by 
2 (l- Vl-Mbo 2 ) is consistent with assumption (c) and leads directly to 
the familiar expression for the Karman-Tsien rule 


°P 


V 

V l-Mco 2 + (l- n/i-Mco 2 ) (cp./2 ) 


(36) 


This rule, in common with the Prandtl-Glauert rule and the present result 
given in equation (28), is termed a pressure correction formula because Cp 
is given in terms of Cp_ and with no further dependence on airfoil 

shape. A comparison of the variations of Cp with indicated by these 
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three relations is shown in 
sketch (f) . A great many- 
other pres sure- correction 
formulas having widely vary- 
ing properties have also "been 
proposed in recent years. 

One that yields results in 
closer agreement with equa- 
tion (28) than the Prandtl- 
Glauert or the Karman-Tsien 
rule has been given by Garrick 
and Kaplan in reference 17. 

A curve illustrating their 
results is included on 
sketch (f) . 


A second important 
method that has been used 
to obtain higher approxima- 
tions for subsonic pressure 
distributions on thin air- 
foils is the method of suc- 
cessive approximation in 
which the solution is 
expressed in a power series in thickness ratio. In this method, the first 
term is the result given by linear theory, and the coefficients of succes- 
sive terms are determined by iteration. Higher approximations cannot be 
expressed in terms of CL. and ML in such a simple and universal manner 

■^i 


Sketch (f) 


as with the pres sure -correct ion formulas, but depend on the airfoil shape 
in a more explicit manner. Although a relatively simple and general pro- 
cedure for the evaluation of the second approximation has recently been 
given by Van Dyke (ref. 18), the determination of the third approximation 
has been accomplished for only a few special shapes. One of these is the 
nonlifting symmetrical circular-arc section for which the second approxi- 
mation has been given by Hantzsche and Wendt (ref. 19) and the third 
approximation by As aka (refs. 20, 21, and 22). 3 Sketch (g) shows a com- 
parison of the variations of Cp with M*, at the midpoint of such an air- 
foil having a thickness ratio t of 0.10, as indicated by equation (28) 
and by the first, second, and third approximations. It can be seen that 
the results obtained by use of equation (28) are identical to those given 
by the first approximation (or linear theory) for small Mach numbers, but 
depart therefrom with increasing Mach number and are much closer to the 
higher approximations for Mach numbers near the critical. It should be 
noted that the curves labeled first, second, and third order represent 
the results indicated by successive approximations to the solution of the 


3 The results for the third approximation given in the present report 
differ from those obtainable directly from the expressions given in either 
reference 20 or 21 and 22 due to the correction of seme misprints. These 
corrections have been verified by correspondence with the author. 
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exact equation for inviscid compress - 
ible flow. It is also of interest to 
compare the results indicated by 
equation (28) with those indicated by 
successive approximations to the 
solution of the simplified equations 
of transonic small disturbance theory. 
The latter results can be calculated 
by use of the following expression, 
which is readily derived from Asaka' s 
result by taking the limiting form 
consistent with the approximations of 
transonic flow theory: 



8 t Ao is (r+i-)Ho g a 

* ^ " 2 ' 0--O 2 T 



f ln2 ) + Kf Ui, - 9 ° 8 ) 


( 7+1) 2 Mqq 4 
(1- J ^, 2 ) r/2 


= - 2.5465 — — — 


0.5132 


( 7 + l)Hx, 2 a 

(l-^o 2 ) 2 T 


0.6339 


(y+l) 2 ^ 4 a 

(i-H=o 2 ) 7/2 T 


(37) 


It can be seen that the curve in 
sketch (h) indicated by equation (28) 
Is somewhat higher than even that 
representing the third approximation, 
but evaluation of its accuracy 
remains difficult because neither the 
exact solution nor an upper bound for 
the results is provided by the clas- 
sical method of successive approxima- 
tion. Attention is called to the 
fact that recent developments in 
transonic flow theory permit the 
establishment of an upper bound by 
application of an alternative method 
of successive approximation that 



Sketch (h) 
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involves the solution of quadratic, rather than linear, equations at each 
step of the iteration process. This process, "based on the methods employed 
in reference 9 , is described in the Appendix, and additional results are «- 

given for the specific case considered in sketch (h) . 


PLOWS WITH FREE -SCREAM MACH NUMBER NEAR 1 


The analyses of supersonic and subsonic flows given in the preceding 
sections have started by introduction of a symbol A for the coefficient 
of cp^ and the assumption that A is nonsingular and varies suffi- 
ciently slowly that it can be regarded as a constant in the initial stages 
of the analysis. Since the results so obtained terminate if A=0, or 
physically if sonic- velocity occurs in the flow field, it is immediately 
clear that some change is necessary to study flows with free-stream Mach 
number near 1 where the transition from subsonic to supersonic flow is an 
essential feature . The technique adopted is to introduce the symbol Ap 
as an abbreviation for the coefficient of <P X rather than thus 


>P = Hxi 2 —■ k ^ (38) 

^ Six 


whence equation (3) may be written as follows: 


<p zz - Vp x » -(i-^ 2 )^ - fp (39) 


If attention is confined to flows with free-stream Mach number 1 so 
that the right-hand side of equation (39) vanishes and if Ap is replaced 
by a constant, the resulting relation given by equation (39) is a linear 
partial differential equation of parabolic type that is familiar from the 
study of one-dimensional unsteady conduction of heat. If approximate 
solutions for flows with free-stream Mach number 1 are sought in this way, 
the analysis proceeds through considerations that are generally applied 
to parabolic differential equations, and is, in some sense, intermediate 
between the mixed elliptic-hyperbolic type of the transonic equation. 

The idea of using the equation for heat conduction for the study of tran- 
sonic flows in this way is not new, but apparently originated with 
Oswatitsch, who suggested it to Behrbohm for the analysis of internal 
flows of nozzles (refs. 23 and 24). The same idea has been applied more 
recently by Oswatitsch and Keune (refs. 25, 2 6, and 27) to calculate the * 
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flow around the forepart of slender todies of revolution at free-stream 
Mach number 1, and they have shown that the results are in remarkable 
agreement with those measured on the front half of a circular-arc body 
of revolution. Although the parameter Ap is regarded throughout as a 
constant, and various means are proposed for the selection of an appropri- 
ate value, it develops that the numerical result for the pressures on the 
forward part of typical smooth bodies of revolution depends so slightly 
on the actual choice that almost any reasonable value can be used for Ap. 

If flows with free-stream Mach number different from unity are con- 
sidered so that the right-hand side of equation (39) remains, and if Ap 
is again replaced by a constant, the resulting differential equation is 
linear and is of e lli ptic or hyperbolic type depending on whether the 
free-stream Mach number, rather than the local Mach number, is less than 
or greater than unity. Maeder and Thommen (ref. 28) have suggested that 
this linearized equation, or its counterpart in three dime nsions, be 
applied, together with a new and arbitrary rule for the selection of a 
value for Ap, to calculate the pressure distribution on complete bodies 
of revolution and on airfoils in two-dimensional flows. The selection 
of an appropriate value for Ap is much more critical for these problems 
than for those discussed originally by Oswatitsch and Keune, however, and 
the replacement of Ap by a constant results not only in serious loss of 
accuracy in many applications, but also in loss of certain essential gen- 
eral features of the solution. In general, results obtained by replac- 
ing Ap, or ciu/cix, with a constant appear to be remarkably accurate if the 
resulting values calculated for cJu/8x are, indeed, nearly constant over 
most of the chord. If, on the other hand, chi/cht varies substantially 
over the chord, no choice of a single value for Ap will suffice to pro- 
vide a useful result. This point is developed further in the course of 
the following discussion and in the Appendix. Some criticisms of the 
above procedure, although principally from a different point of view, 
have appeared in a note by Miles (ref. 29). 

In the present analysis, it is assumed once again that Ap is non- 
singular and that it varies sufficiently slowly that it can be considered 
as a constant in the initial stages of the analysis in which a nonlinear 
ordinary differential equation is established for u on the airfoil sur- 
face. The final result for u is determined by integration of this dif- 
ferential equation and restores, to a large degree, the effects of the 
variation of Ap along the chord. The result for the first stage of the 
analysis requires the solution of equation (39) subject to the boundary 
conditions given in equation (l) and can be obtained by application of 
standard procedures. The solution has two distinct forms depending on 
the sign of Ap. The form associated with positive Ap is appropriate 
for application to regions where the flow is accelerating, whereas that 
associated with negative Ap may be appropriate for application to 
regions where the flow is decelerating. The analysis of accelerating flow 
will be developed first. The direct problem in which the airfoil shape 
is specified and the pressure distribution is sought, and the inverse 
problem in which the pressure distribution is specified and the associated 
airfoil shape is sought are discussed for each case. 
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Accelerating Flows , Direct Problems 


Approximate solution of equations for arbitrary airfoil shape .- A 
relation for u at the airfoil surface derived by consideration of equa- 
tion (39) with positive Ap, the boundary conditions stated in equa- 
tion (l), and the form of Green’s theorem associated with the left side 
of equation ( 39 ) (see Appendix) is 



The two alternative expressions for Up are completely equivalent. The 
first is more concise and will be used in the following equations, but 
the second is often somewhat simpler to evaluate. If the free -stream 
Mach number is unity, the double integral vanishes and up can be cal - 
culated directly. The result so obtained corresponds to that found by 
application of Maeder and Thommen' s proposal of reference 28. (it should 
be noted, however, that the general expression for <p given in refer- 
ence 28 is incorrect owing to improper treatment of plus and minus signs.) 

If the free- stream Mach number is not unity, equation (40) is an 
integral equation, and it might appear that little progress toward a 
solution has been made. If attention is confined to the vicinity of the 
airfoil and to Mach numbers near unity, however, it is only necessary to 
approximate cpgg well locally and it is sufficient to substitute Ap/k 
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for or du/B| in tlie double integral. The integral can then he 

evaluated and the following relationship results: 


u = (1 -Hx, 2 )Uqo _ Uco _d r X d| 

p ^(r+i) 

If, once again, kux is restored in place of Xp so that, in effect, the 
local -value is used at each point, and the subscript P on u is dropped, 
a simple nonlinear ordinary differential equation Is obtained for u 



_ (l-y^Uoo 

(r +1 ) 


u„ 


^3tk(du/dx) 


- r 

** Jo 


dZ/d| 

*7 x-i 


d| 


(42) 


Equation (42) can be written in the following form, upon rearranging 
terms and squaring both sides. 


u 


(l-Qu,. I 2 du = Uco 3 

Moo 2 (7+1) J dx nMoo 2 (r+l) 



the 


(43) 


As in the other cases, equation (43) can be solved readily by separation 
of variables, and the constant of integration can be evaluated by intro- 
duction of the additional condition that equation (4l) provides the cor- 
rect location x=x* for the sonic point, or point where 
u=(l-Moa 2 ) Uoc/Moo 2 ( 7+1) . The result is 


L u 


(i-m 00 2 )u 00 -| 3 
^=o 2 (7+1) - 



(44) 


where x* is the value for x for which 


d 

dx 



dZ/d| 

7^7 


d| = 0 


(45) 


This method of evaluation of the constant of integration Is completely 
analogous to that employed in the analysis of subsonic and supersonic 
flows, and Is necessary in the analysis of flows with free-stream Mach 
number near 1 in order to avoid infinite pressure gradients at the point 
of sonic velocity on smooth airfoils. This method, moreover, provides a 
mechanism for the introduction of direct upstream dependence on airfoil 
shape in the subsonic region, and its exclusion in the supersonic region. 
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The corresponding relation for the pressure coefficient Cp is obtained 
by combination of equations (2) and (44) and is 


, -2(l-^o 2 ) 2 3 

P H>o 2 (7+1) mV (7+1) 



(46) 


An alternative expression in terms of the transonic similarity parameter 
is 


- _ [V( 7 +l )] 1/3 

^ T 2 / 3 


Gp= 2| ro -2 
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d(z/r)/dj 

•J 


V V 3 




dxi 


= 2g + 

3 oo 


'£ «0 
=00 


(47) 


The variation of Cp with expressed by equation (47) is exact, within 

the approximation of transonic small disturbance theory, for flows with 
free-stream Mach numbers very near unity, and is associated with the fact 
that the local Mach number distribution on an airfoil is independent of 
the free-stream Mach number at values of the latter near unity. This 
phenomenon has been discussed previously in references 30 , 31 * 6, and 
elsewhere . 


Once Cp is known, the pressure drag d can be readily calculated 
by use of the following relation 


c d s 


P a 
-f Uco C 



(48) 


Application to single-wedge airfoils .- Sufficient theoretical and 
experimented, results are available at the present time to provide consid- 
erable insight into the accuracy and usefulness of equation (46) or (47) • 
The shape for which the greatest amount of information is available is 
the single-wedge profile for which both theoretical and experimental pres- 
sure distributions are available. Thus, con- 
sider a single-wedge profile of maximum 
thickness t and chord c/2 as illustrated 
in sketch (i) . The ordinates of the airfoil 
upper surface are 

Sketch (i) 
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* 


Z = t — » tx for 0 <x< c/2 

c 

Z = t/2 = tc/2 for x> c/2 




and the semiapex angle d is equal, to the order of accuracy of thin, air- 
foil theory, to r. Substitution of equation (49) for Z into equa- 
tions (45) and (46) provides that the sonic point is at the shoulder 
(x*=c/ 2 ) and that the pressure distribution on the surface of the wedge 
at free-stream Mach numbers near 1 is 


IX* 2 ( 7 + 1 )] 1/3 

C p S t 2/3 °p 


26oo-2 




(50) 


A plot of the results for Mach 
number 1 is shown in sketch (j) 
together with the corresponding 
theoretical results given by 
Guderley and Yoshihara in refer- 
ence 32. Although some approxima- 
tions are introduced in the course 
of the latter analysis, the results 
are generally regarded as virtually 
an exact solution of the equations 
of transonic small disturbance 
theory. Also included in sketch ( j) 
are experimental results for M^=l 
obtained in the Langley annular tran- 
sonic wind tunnel and reported by 
Habel, Henderson, and Miller in ref- 
erence 33 * 

Since the comparisons shown in 
sketch ( j) indicate that equa- 
tion ( 50 ) provides an approximate 
solution for the pressure distribu- 
tion on a single-wedge profile 
at ^*=1 that is probably satisfac- 
tory for most purposes, and since 
the variation of Up with | TO given 
by equation (47) or ( 50 ) is exact, 
within the .framework of transonic 
small disturbance theory, at 1 ^= 1 , 
the principal question remaining in 
the evaluation of the degree of 
approximation afforded by use of 
equation ( 50 ) is to define the range 
of 1 ^, or Mach number, over which 
it applies. Accordingly, sketch (k) 
has been prepared to summarize the 



Sketch (k) 
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results given previously by Guderley and Yoshihara for Mach number 1 
(ref. 32), by Vineenti and Wagoner for Mach numbers slightly greater than 
unity (ref. 34), and by Yoshihara for Mach numbers slightly less than 
unity (ref. 35) • The latter two sets of results were obtained by lengthy 
numerical calculations and, together with the results for Mach number 1, 
are generally regarded as being very close to those that would be given 
by exact solutions of the equations of transonic small disturbance theory. 
The results are plotted in terms of Cp-2^ so that the pressure distri- 
butions for Mach numbers very near unity should determine a single line. 

It can be seen that the variation of "Sp with g^ indicated by equa- 
tion (47) or (50) holds until the absolute value of g ra is nearly one- 
half. At greater values, the results begin to tend toward those asso- 
ciated with purely subsonic or purely supersonic flows, and equations (47) 
and (50) are no longer applicable. 

Experimental measurements of the flow around single-wedge profiles 
at free-stream Mach numbers both less than and greater than unity have 
been made by Liepmann and Bryson and reported in references 30 and 31- 
Results were obtained for three different profiles having semiapex angles 
of 4-1/2°, 7-l/2°, and 10°. Plots of the experimental pressure distri- 
butions for the test Mach numbers closest to unity for each profile are 
shown in sketch (z) together with the theoretical pressure distribution 



Sketch (z) 

calculated by use of equation (50). Additional experimental data for 
other Mach numbers are not included on sketch (2) since those shown are 
already for values of g^ that are somewhat outside the range of valid- 
ity of equation (50). Only the theoretical results given by equation (50) 
axe included since examination of the data reveals that these results dif- 
fer less from the theoretical results shown on sketch (k) , for Mach num- 
bers near unity, than the experimental curves differ from either set of 
theoretical curves, or even than the experimental curves differ among 
themselves. Perhaps the most prominent discrepancy is that which occurs 
near the shoulder. Theory Indicates that sonic velocity (^p-2g co =0) 
occurs at the shoulder, whereas the experimental data, particularly that 
of Liepmann and Bryson, consistently indicate that sonic velocity occurs 
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forward of the shoulder. It is interesting to observe that this 
discrepancy is greatest for the thinnest airfoil tested. 


The foregoing results may be contrasted with those obtained by 
direct use of eq uatio n (4l) in which case Cp- 21 ^ is found to be pro- 
portional to iA/VL It is clear from this comparison that there is no 
constant with which Ap can be replaced that would provide a satisfac- 
tory solution for the pressure distribution on a single-wedge airfoil at 
free-stream Mach numbers near unity. 

The pressure drag of single-wedge profiles at Mach numbers near 1 
can be found easily by integration of equation (48) with the relations 
given by equations (49) and (50) substituted for Z and Cp. The result 
is 


[^(7+1) ] 1/S e d - 

= s Cfl = 


t 5/3 


G d = 2 | m + 2 


2 ) 


(* 


1/3 


(!) - 2 5oo + 1-758 (51) 


where r represents the gamma function. This result compares very 
favorably with the value of 1-75 fox* or Mach number 1, given by 

Cole in reference 36 as that obtained by numerical integration of the 
pressure distribution given by Guderley and Yoshihara in reference 32. 
Cole’s own theory for the drag of a single-wedge airfoil at high subsonic 
speeds , which is fundamentally somewhat less accurate than Guderley and 
Yoshihara’ s theory, gives the value I .67 for 1^=0. 

Application to biconvex circular-arc airfoils .- In addition to data 
for three single-wedge profiles, Bryson also gives in reference 30 exper- 
imental results for transonic flow around the front half of an 8 . 8 -per cent 
thick biconvex circular-arc airfoil followed 
by a straight section an illustrated in 
sketch (m) . Since the pressure distribution 
on the curved portion of this profile is the 
same at Mach numbers near unity as that on 
the front half of a complete circular-arc 
airfoil having the same thickness ratio, and 
additional experl mental data are available 
for the latter airfoil although for other 
thickness ratios, the following analysis is 

developed for a complete biconvex circular-arc airfoil. It is moreover 
sufficient, in thin airfoil theory, to approximate the ordinates Z for 
a circular-arc airfoil by those for a parabolic-arc airfoil, thus 

t = ^ [f - CsT ] 



( 52 ) 
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z=2tc[£-<£> 2 ] 

r 


W/mna-y 


Sketch (n) 


where r is the thickness ratio as 
indicated in sketch (n) . Substitution 
of this relation for Z into equa- 
tions (4-5) and (-4-7) yields the following 
result for the pressure distribution on 
the airfoil surface. 


77 _ [M^ 7+1) 1 1/3 

°p 


Cp = 21^-2 


f K 4 f) 




(53) 



Sketch (o) 


Sketch (o) shows a comparison 
of the pressure distributions for 
Mach numbers near unity calculated 
by use of equation (53) with those 
obtained from Bryson's experiments 
with the half airfoil. As for the 
■c single-wedge profile, the results 
are again plotted in terms of 
Cp-21^ since experimental results 
are available only for Mach numbers 
somewhat different from unity. It 
can be seen by comparison of 
sketches (l) and (o) that the theo- 
retical and experimental results are 
in much better agreement for the 
front half of a circular-arc airfoil 
than for single-wedge profiles. 
Experimental pressure distributions 
for transonic flow past four complete 
biconvex circular-arc airfoils having 
thickness ratios of 6, 8, 10, and 
12 percent have been given by Michel, 
Marchaud, and Le Gallo in refer- 
ence 37. Their results for Mach 
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number 1 are plotted In 
sketch (p) together with 
the theoretical results 
calculated by use of equa- 
tion (53) . These results 
are presented in terms 
of Gp because transonic 
theory indicates that the 
pressure distributions for ^ 
all four airfoils should 
then define a single curve 
independent of the thick- 
ness ratio. Results for 
Mach numbers other than 
unity are not included on 
this plot because the 
variation of TJp with 
for small g„ indicated 
by equation (47) and sub- 
sequent relations is not 
only simple but is anply 
verified by the preceding 
comparisons and by similar 
discussions elsewhere in the literature of transonic flow. It can be 
seen that the theoretical and experimental results are in substantial 
agreement. The most notable discrepancy is that found near the trailing 
edge, and can be attributed to flow separation induced by boundary- layer 
shock-wave interaction. It can also be seen that the agreement between 
theory and experiment is not so good for the complete airfoils, particu- 
larly the thinner ones, as for the half airfoil. Part of the discrepancy 
for the complete airfoils may possibly be attributed to the experimental 
technique in which the airfoil is simulated by a bump on the tunnel wall 
and is hence imbedded in the wall boundary layer. Some comments on this 
method of testing have been given recently by Carroll and Anderson in 
reference 38. 



The pressure drag of circular-arc airfoils at Mach numbers near 1 
has been found by numerical integration of equation (48) with the rela- 
tions given by equations (52) and (53) substituted for Z and Cp and is 


tn/(ri)] v * a ,- a , ii.7 7 

T 5/3 


(54) 


The pressure drag of the front half can be evaluated in the same way by 
changing the upper limit in equation (48) to c/2 and is 


c d = 2£ ro + 1.13 


(55) 


The integrations required to determine the drag results given in 
equations (54) and (55) were evaluated numerically using Simpson’s rule 
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together with an analytic determination of the contribution of the region 
in the immediate vicinity of the leading-edge singularity. Sufficiently 
fine intervals were used that the resulting values are estimated to be 
accurate to within about one digit in the third significant figure, as 
judged by comparison with the results of similar calculations made with 
wider intervals. It is necessary to use very fine intervals, particu- 
larly near the nose , to achieve such accuracy, and intervals as small 
as 0.00005c were used in some cases. 

Application to a family of airfoils having the point of maximum 
thickness displaced aft of the midchord station .- The primary object of 
this present section is to present some comparisons of calculated and 
measured pressure distributions at Mach number 1 on a n umb er of specific 
airfoils that have the point of maximum thickness aft of the midchord 
station. The experimental data are from reference 39 by Michel, Marchaud, 
and Le Gallo, and are for members of the family of airfoils having ordi- 
nates given by 


f = A [f-Cl) n ] <*> 

where A and n are constants for each airfoil and n is greater than 
unity. The values selected for A and n determine the thickness 

ratio t and the location (x/c)„ of the point of maximum thickness 

“max 

according to the relations 


(57) 


= (58) 
n n/(n-i) 


The biconvex circular- arc airfoils discussed in the preceding section are 
special cases of the present family that correspond to n=2. The point 
of maximum thickness is located forward or aft of the midchord station 
depending on whether n is less than or greater than 2. The particular 
airfoils tested by Michel, Marchaud, and Le Gallo are special cases that 
correspond to either n=3.38 or 6.05 and have the point of maximum thick- 
ness at 0.60 or 0.70 chord. As in the earlier work by the same investiga- 
tors on biconvex circular-arc airfoils, the results were obtained by 
simulating the airfoil by a bump on the tunnel wall and are again subject 
to criticism regarding the influence of the wall boundary layer. 
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Substitution of equation ( 56 ) for Z into equations (45) and (46) 
yields the following result for the pressure distribution on the airfoil 
surface for Mach numbers near unity: 



[H/( 7+1)1 1/3 

T2/3 




= 21^-2 


_3_ 
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n n/(n-i) 


>1/8 
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(n-l)r(n-i) 


. /• v 11-1 

m 


a/s 


ac[r(n+i) 3 


2 / \2(n-i) 
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x \ ~f\ c / n -1 


2 (n-l) r(n-|^ 


in J* . r fa +1 I + --3- 


, 1 A 2 (n-l) 

r( n " g 


(59) 


where r represents the gamma function. If m is any positive integer 
greater than unity, the following relations are useful for evaluation of 
the gamma function: 

f(m) = 1 - 2-3 - - - (m- 1 ) = (m-l)I , r(l/2) = *«/* 1 



Substitution of n=2 in equation (59) reproduces equation ( 53 ) for 
the pressure distribution on biconvex circular-arc airfoils as a special 
case. Equation ( 50 ) for the pressure distribution on a single-wedge air- 
foil can also be obtained as a limiting case by setting n=<» and noting 
that the chord of the wedge is designated here as c rather than c/2. 
Theoretical pressure distributions on the airfoils tested by Michel, 
Marchaud, and Le Gallo are obtained from equation (59) by substitution 
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Sketch (r) 


of the values 3.38 and 6.05 
for n. Sketches (q) and 
(r) show comparisons of the 
theoretical and experimental 
pressure distributions for 
1^=1 for the two groups of 
airfoils. Except for the 
discrepancy near the trail- 
ing edge which can again be 
attributed to boundary- 
layer shock-wave inter- 
action, it can be seen that 
the theoretical and experi- 
mental results are in at 
least qualitative agree- 
ment. Some differences 
occur, however, in the 
levels of the pressure dis- 
tribution curves. Compari- 
son with the results for 
the circular-arc airfoils 
shown in sketch (p) reveals 
that the same trend is in 
evidence for those air- 
foils, although to a lesser 
degree, and that the differ- 
ence between the theoretical 
and experimental results 
increases as the point of 
maximum thickness moves 
rearward. It is not clear 
at the present time whether 
this discrepancy is to be 
attributed principally to 
the shortcomings of the 
theoretical or the experi- 
mental results. 

Application to a 
family of airfoils having 
the point of maximum thick - 
ness forward of midchord 
station .- The test program 
of Michel, Marchaud, and 
Le Gallo reported in ref- 
erence 39 and discussed in 
the preceding section also 
included tests of each of 
the airfoils reversed in 
the wind tunnel so that the 
point of maximum thickness 
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is located forward of the midchord station. The particular airfoils 
tested are thus specific cases of the family of profiles described by 
equation ( 5 6 ) with x/c replaced by l-(x/c) ; that is, the ordinates are 
given by 


Z 

c 


A 



x 

c 



(61) 


where A and n are again constants for each airfoil and n is greater 
than unity. The values selected for A and n determine the thickness 
ratio r and the location ( x / c )z max P°4Q.t of maximum thickness 

according to the relations 


1 




( 62 ) 


T - 2A(n-l) 
n n/(n-i) 


(63) 


Biconvex circular-arc airfoils are special cases of the present family 
that correspond to n=2. The point of maximum thickness is located forward 
or aft of the midchord station depending on whether n is greater or less 
than 2. The particular airfoils tested by Michel, Marchaud, and Le Gallo 
are special cases that correspond to either n=3 . 38 or 6.05 and have the 
point of maximum thickness at 0.40 or 0.30 chord. 

Since the integrations encountered when equation ( 6 l) is substituted 
in equations (45) and (46) for the determination of Cp are more diffi- 
cult than those encountered in any of the preceding exangles, no general 
formula will be given for arbitrary n. Two formulas of more restricted 
generality are given, however. One is applicable when n is any integer 
greater than unity, the other when n is one-half plus any positive 
integer. The first iB 
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( 64 ) 


■where the symbol f x# is used to denote the difference between the pre- 
ceding expression with first x and then x* substituted for x. that 
is. 


f ( x ) = f(x)-f(x*) 

where x*, the location of the sonic point, is found from 

n-i 

n-l+r(n+l) o/jr ^ (-x*/c) _ Q 

v=i r(n-v)T (v +i ^ 

and T refers again to the gamma function. The second is 
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where B = 




/ 1> 


, and F(£) = 
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Attention is called to the fact that v and p are positive integers so 
that when n=(3/2) all the summations drop out. 


Again x* is the location of the sonic point and it is found from 





Substitution of n=2 
in equation (64) again 
reproduces equation ( 53 ) 
for the pressure distribu- 
tion on biconvex circular- 
arc airfoils. Sketches (s) 
and (t) show comparisons 
of the pressure distribu- 
tions measured at Mach num- 
ber 1 by Michel, Marchaud, 
and Le Gallo with those 
calculated by use of equa- 
tions (64) and ( 65 ) . The 
experimental results shown 
in sketch (s) are for air- 
foils that have the point 
of maxi mum thickness 
located at 0.40 chord cor- 
responding to a value 
for n of 3 . 38 . Since 
results could not be cal- 
culated analytically for 
this value for n, theo- 
retical results are shown 
for both n=3.0 and n=3 . 5 • 
The corresponding locations 
for the point of maximum 
thickness can be readily 
calculated using equa- 
tion ( 62 ) and are 0.423 
and 0.394 chord. Similarly, 
the experimental results 
shown in sketch (t) are 
for an airfoil that has 
the point of maximum thick- 
ness located at 0.30 chord 
corresponding to a value 
for n of 6 . 05 ; whereas 
the theoretical results 
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are for airfoils tliat have the point of maximum thickness located at 
O. 3 OI chord corresponding to a value for n of 6.00. These results con- 
tinue the trend noted in the preceding section that the agreement "between 
the present experimental and theoretical results improves as the location 
of the point of maximum thickness moves forward along the chord. The 
principal discrepancy remaining is , in fact, reduced to that near the 
trailing edge associated with boundary- layer shock-wave interaction, and 
is therefore beyond the scope of any inviscid theory. 


Accelerating Flows, Inverse Problems 


Approximate solution of equations for given pressure distribution . - 
Although all of the preceding discussion is concerned with the calculation 
of the pressure distribution on an airfoil of specified geometry, an 
equally important problem in many engineering situations is the design of 
an airfoil to have a specified pressure distribution at a given Mach num- 
ber. This poses no new problem in the analysis of purely subsonic or 
purely supersonic flows by the present methods because the specification 
of Cp permits the determination of Cp^ through application of equa- 
tion ( 31 ) or ( 17 ) , and the inversion problem is reduced to the familiar 
inversion problem of linear theory. The necessary relation for flows with 
free-stream Mach number near 1 can be derived from consideration of equa- 
tion (42) as an integral equation in which u and du/dx are given and 
the unknown appears in the integrand of a definite integral. This equa- 
tion can be inverted readily since it has the same form as the relation 
encountered in the solution of Abel’s integral equation (see, e.g., 
ref. 40, pp. 483-484). The inversion thus has the form of Abel’s integral 
equation, and is the following in the present application: 


dz = _ _i /k r x f u 1 /~&u dg 

dx Uoo V it J Q L ^ 2 (7+ 1 ) -I* J x— £ 


( 66 ) 


The desired relation for the ordinates Z(x) of the airfoil can now be 
found by a second integration, and is the following if it is assumed 
that Z is zero at the leading edge (x = 0) . 


Z = Z* - 
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where Z* is the ordinate at the point x* where u is zero. It is 
interesting to note that the two alternative expressions for Z lead 
to identical results although the apparent regions of dependence, as 
indicated by the limits on the integrals, are quite different. The same 
result expressed in terms of Cp or Cp is 


Z = 


1 1 ^(7+ 1) r x 

2 si 2* J 0 


dxi 






dg 


( 68 ) 


or 


Z = 







(69) 


A simple application or check of these relations is furnished by substi- 
tution for Cp of the relation given in equation (50) for single-wedge 
profiles, whence Z is found to be equal to 9x between x=0 and x=c/2. 
In the same way, substitution of equation (53) for Cp leads directly 
to equation (52) for the ordinates of a circular-arc airfoil, etc. 

Application to airfoils with constant pressure gradient .- An example 
that permits an additional comparison with an, existing theoretical result 
given by Guderley in reference 4l, is furnished by consideration of the 
problem of determining the shape of an airfoil having a constant negative 
pressure gradient at Mbo=l, thus 

= = <*» 

where A = -dCp/d(x/c) is a positive constant. Substitution of this 
relation for Cp in equation (68) and integration leads directly to the 
following result: 


Z = 2 / 2 +I A s/a ( x \ 3/2 ( x* _ 2 x \ 

0 3 v 2jt V c / V c 5 c / 


(71) 


The special case considered by Guderley is obtained by insertion of the 
values x*/c = 3/4 and A = 6/5, whence equation (71) reduces to 
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Sketch (u) shows a scale drawing of 
the profile calculated by use of 
equation (72) and of that given by 
Guderley. It is evident from the 
sketch that the present theory indi- 
cates larger values for the ordi- 
nates Z than are given by Guderley. 
Although the latter results are given 
only in graphical form, and are hence 
difficult to determine with precision 
the two sets of values for Z appear 
to be related by a constant ratio of 
approximately 9 to 8. 



Sketch (u) 


The case considered by Guderley and discussed above results in a 
shape that does not close at the stem. It can be seen immediately from 
equation (71) , however, that a closed airfoil will result if x*/c is 
equated to 2/5, in which case equation (71) reduces to 


Z 
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Jt M a 3/2 

15 v 2* 


or (-0 


(73) 


A plot of the results is shown in 
sketch (v) . 



Decelerating Plows 


Although it is clear that the preceding relations are not appropriate 
for the analysis of flows with free-stream Mach number near 1 that decel- 
erate smoothly through sonic velocity, it might appear that the proper 
expressions could be derived by formal application of the procedures des- 
cribed in the preceding section for positive Ap to the approximate solu- 
tion of equation (39) for negative Ap. The analysis leads, in the absence 
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of contributions from shock waves, to the following relation for u at 
the airfoil surface instead of equation (40) 


U P 


tka f2pu + is r at 

^-irAp ^ o x n/I-x Ap Sx J,*, 



(74) 
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The principal difference between the results for the two cases is that 
the value of u(x) indicated by equation (74) depends on conditions down- 
stream of the point x, whereas that indicated by equation (40) for posi- 
tive Ap depends on conditions upstream of x. This difference is a 
fundamental property of equation ( 39 ) and necessitates a change in the 
argument required to disregard additional contributions from shock waves, 
because now it is the oblique shock waves situated downstream, rather than 
upstream of x, that furnish a contribution to u(x) . If, however, there 
are no oblique shock waves downstream of x, or if the contributions 
resulting from additional Integrals over the shock surfaces are disre- 
garded so that equation (74) can be used as a starting point in the anal- 
ysis, the following result is obtained by proceeding in a fashion 
analogous to that employed in the derivation of equation (47) from (40) : 
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The symbol x* again refers to the location of the sonic point and is 
equal to the value for x for which 


-i rsastdi-o 


(76) 


The corresponding relation for the shape of an airfoil associated with a 
given pressure distribution at some free- stream Mach number near unity 
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can also be found and is the following, again assuming that Z is zero 
at the leading edge: 


Z 




(Cp-20 Ja-CpM 




(77) 


No further use is made in this paper of equations (74) through (77) 
for decelerating flows at free-stream. Mach numbers near 1. As will become 
more evident at a later point in the present discussion, it would appear 
necessary to use such formulas for the analysis of flows decelerating 
through sonic velocity, but the region of dependence in these relations 
is such as either to cast suspicions on their applicability or to require 
the occurrence of exceptional coincidences. On the other hand, two- 
dimensional flows that decelerate smoothly through sonic velocity appear 
to be very exceptional physically. Further investigation is needed before 
additional remarks can be made regarding the role of the parabolic case 
with negative Ap in the analysis of transonic flows. 

The next section will be concerned with an alternative analysis of 
certain cases in which decelerating flows occur. 


Combination of Accelerating and Decelerating Flows 


The calculation of the pressure distribution at Mach numbers near 1 
on an airfoil having such a shape that the velocity increases over part 
of the chord and decreases over the remainder cannot be accomplished by 
direct application of any of the relations developed in the preceding 
sections. On the one hand, the parabolic method described for flows with 
free-stream Mach number near 1 permits the analysis of flows that pass 
through sonic velocity, but fails when the velocity gradient is zero. On 
the other hand, the elliptic and hyperbolic methods described for subsonic 
and supersonic flows permit the analysis of flows with zero velocity gra- 
dient, but fail if the local velocity is sonic. The breakdown in each 
case is associated with the fact that the basic partial differential equa- 
tion for each case, that is, equation (8), (21), or (39 ) , assumes a degen- 
erate form when A iB zero. Such cases are, nevertheless, interesting 
and important since they can occur in practical applications, and the pre- 
sent section is concerned with their discussion. The procedure adopted 
is based on the idea of joining together various of the results derived 
in the preceding sections in such a way that the failings associated with 
vanishing A are avoided, rather than on a complete re-analysis of the 
problem from a sufficiently general point of view to encompass the entire 
problem in a single sweep. 


In order to fix the ideas, consider the problem of calculating the 
pressure distribution at Mach number 1 on the airfoil with cusped trailing 
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edge illustrated in sketch (w) for 
which experimental data are available 
from reference 42 by Michel, Marchaud, 
and Le Gallo. The front half of this 
airfoil is the same as that of a 
biconvex circular-arc airfoil having 
a thickness ratio of 0.10, but the rear 
half is shaped in such a manner that 
an inflection point is located at 
0.75 chord and that the trailing-edge 
angle is zero. The ordinates and 
slopes of the rear half of this air- 
foil are shown graphically in 
sketch (x) . No analytic egression 
is given in reference 42 for the ordi- 
nates of the rear half. 


The pressure distribution on the 
front part of the airfoil where the 
passage through sonic velocity occurs 
can be calculated by use of equa- 
tions (45) and (46), since it does not 
depend on the shape of the rear half of the airfoil. This means that the 
pressure distribution on the front half of the airfoil described above is 
given specifically by equation (53) for x/c between 0 and 1/2. It Is 
clear that the pressure distribution on the entire rear half of the air- 
foil cannot be determined by use of equations (45) and (46) because the 
results so calculated indicate a point of zero pressure gradient in the 
vicinity of the inflection point. Although this detail, in itself, is 
not incorrect, it signals the breakdown of the parabolic method that 
occurs when A is zero. Positive evidence of the breakdown is provided 
by the fact that the calculated pressures decrease downstream of the point 
of zero pressure gradient rather than increase as indicated by the experi- 
mental data shown on sketch (w) or by simple considerations of supersonic 
flow. These results, furthermore, cannot be joined to those obtained by 
use of equations (75) and (l6) for the part of the airfoil downstream of 
the point of zero pressure gradient because the two sets of equations do 
not indicate the same location for this point. This situation should not 
be too surprising since the procedures should not be expected to fail 
abruptly when A is precisely zero, but gradually as A approaches zero. 

There exists another possibility for the determination of the pres- 
sure distribution on the rear half of the present airfoil by joining 
together solutions. It Is to use the formulas developed for supersonic 
flow, but with the final constant of integration adjusted so that the 
pressure is equal, at the point of connection, to that given by the solu- 
tion for the forward part of the airfoil. This procedure corresponds to 
the use of simple wave theory for the calculation of the difference in 
pressure between an arbitrary point on the rear of the airfoil and the 
point of connection. In this way, the following equation results for 
the pressures on the rear of the airfoil at Mach numbers near unity: 
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where Z* ■ refers to dZ/dx, and Cp(X) is the value of Cp at x=X. 

The most convenient point for joining the two solutions in the appli- 
cation described in sketches (w) and (x) is at x/c = 0.50. Then the 
pressures on the forward half of the airfoil can he calculated directly 
by use of equation (53), the values of Cp(X) and Z’ (X) are 


Cp(X) - Cp ( f ) - 2^-2 | (-1 + in 4) 


1/3 



and the following expression results for the pressures on the rear half 
of the airfoil at Mach numbers near unity: 



L^x) 2 (y+i) ] 1/3 

1-2/3 


°P = 2 ^- 2 



1 + in 4) 


- _X ^ 

2t dx 


1 2/3 


(80) 


Sketch (y) shows a compari- 
son of the experimental 
pressure distribution for 
Mach number 1 given by 
Michel, Marchaud, and 
Le Gallo in reference 42 
and the corresponding theo- 
retical values calculated 
using equations (53) and 
(80) together with the 
values for dZ/dx given 
in sketch (x) . The theo- 
retical and experimental 
results bear about the 
same relationship to each 
other as those shown pre- 
viously for biconvex 
circular-arc airfoils 
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Sketch (y) 


although effects of boundary- layer shock-wave interaction extend over a 
larger fraction of the chord of the cusped airfoil. This difference is 
in agreement with the results obtained from schlieren photographs and 
given in reference 42 that indicate that the shock wave meets the airfoil. 
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at Mach number 1, at 78-percent chord, for the cusped airfoil and at 
95-percent chord for the biconvex circular-arc airfoil of the same 
thickness ratio. 


It is apparent that the pressures computed over the rear half of the 
airfoil by using equation (80) will tend to be somewhat too negative 
because the use of this relation corresponds to the use of simple wave 
theory and hence disregards the influence of a family of incoming compres- 
sion waves arising from the sonic line. Some idea of the magnitude of 
this effect can be gained by examination of sketch (z) which shows a 

comparison of the pressure distribu- 



tion on biconvex circular-arc air- 
foils at Mach number 1 calculated 
using equation (53) for the entire 
airfoil, with those calculated using 
equation (78) for various fractions 
of the chord. 

A further example involving 
accelerating and decelerating flows 
at Mach number 1 is furnished by 
examining the case of the symmetrical 
double-wedge airfoil of arbitrary 
thickness ratio for which a solution 
has been given by Guderley and 
Yoshihara in reference 32. 

Sketch (a 1 ) shows a plot of their 
result together with the correspond- 
ing result calculated by the proce- 
dures described above. The result 
for 0<x<c/2 is calculated by 
use of' equation (50)- That 
for c/2<x<c is calculated by use 
of the following equation which is 
obtained from equation (78) by 
equating Cp(X) to 0, Z 1 (X) to t, 

and Z r (x) to -r. 

Cp = -2(3) 2/3 (81) 


. The difference between the two pres- 

Sketch (a 1 ) sure distributions on the rear half 

of the airfoil is again the result of the neglect, in the present analysis, 
of the contribution of the family of incoming compression waves arising 
from the sonic line. 


It is evident from these three examples that the present procedure 
of joining solutions is capable of yielding results that are at least 
qualitatively correct and that, although somewhat greater accuracy is 
undoubtedly desired, this procedure may be useful until such time as a more 
general analysis of accelerating-decelerating flows is accomplished. 
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COMPARISON OF RESUITS FOR MACH NUMBER 1 WITH 
THOSE FOR OTHER MACH NUMBERS 


The remainder of this report is concerned with the discussion of some 
selected examples that illustrate the relation between results for Mach 
number 1 and those for other Mach numbers. This discussion is divided 
into two parts. The first is concerned with comparisons of pressure dis- 
tributions on the same airfoil at different Mach numbers, and the second 
with pressure drag. 


Pressure Distribution 


Attention is directed in this section to comparisons, for a number 
of airfoils, of the theoretical pressure distribution for Mach number 1 
with that for the highest Mach number for purely subsonic flow and that 
for the lowest Mach number for purely supersonic flow. Pressure distri- 
butions for these two Mach numbers, designated more briefly as the lower 
and upper critical Mach numbers, respectively, are of particular signifi- 
cance not only because they represent the results associated with the 
bounds of the transonic range, but also because they are typical of the 
pressure distributions for all purely subsonic or purely supersonic flowB. 
All of the theoretical results shown are calculated by application of the 
general expressions derived in the present paper. The results for Mach 
number 1 are the same as discussed in detail in the preceding section. 
Subsonic pressure distributions are calculated by use of equation (28) , 
and supersonic pressure distributions 
by use of equation ( 15 ) . The lower 
critical Mach number is determined 
from equation ( 30 ) by replacing Cp -2 


-.3 


L cr 


with the most negative value of 


and 




, solving 
critical Mach 


that occurs in each case 
for M^. The upper 
number is determined from equa- 
tion ( 16 ) by replacing (dZ/dx) cr 
with the value of dZ/dx at the 
leading edge and solving for M^. 

Consider, first, the single-wedge 
airfoil for which the pressure distri- 
bution at Mach number 1 Is given by 
equation ( 50 ) and illustrated graphi- 
cally in sketch (j). Sketch (b’) 
shows a comparison of this result for 
the specific case of a wedge having a 
semiapex angle Q of 0.10 radians 
with those for the upper and lower 
critical Mach numbers. The lower 
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critical Mach number is, of course, zero because the velocity is sonic at 
the comer for all free-stream Mach numbers less than the upper critical. 
The pressure distribution for Mach number 0 is given by 


°p = S - - f 


In 


2x/c 

l-2x/c 


(82) 


-3 


-2 


-I 





a 

a 

■ 



a 

a 


55 

1 

|Qf 





mm 




93 

R 


■H 

gs 

3 

a 

B 

a 


■ 

a 




a 






Sketch (c*) 

pressure distribution for the lower 
of the following expression for Cp 
the auxiliary relation of equation . _ 

sion given in equation (52) for 0<x<c/2 and by rc/2 for x>c/2: 


The three curves shown on sketch (b* ) 
suffice to show that the pressure 
distribution on a single-wedge air- 
foil at Mach number 1 bears a much 
closer resemblance to that at the 
lower critical Mach number than to 
that at the upper critical Mach num- 
ber. It is interesting to note, 
moreover, that the difference between 
the pressure distribution at Mach 
number 1 and that at Mach number 0 
is very nearly constant across the 
chord. 

Consider, next, the half 
circular-arc airfoil for which the 
pressure distribution at Mach num- 
ber 1 is given by equation (53) and 
illustrated graphically In 
sketch (o) . Sketch (c 1 ) shows a 
comparison of this result with those 
for the upper and lower critical Mach 
numbers. The computation of the 
critical Mach number involves the use 
which Is obtained by integration of 
(28) with Z replaced by the expres- 






2x/c " 
l-2x/c 


(83) 


The results illustrated in sketch (c 1 ) display a remarkable property that 
the subsonic part of the pressure distribution at Mach number 1 differs 
from the pressure distribution at the lower critical Mach number by nearly 
a constant, and that the supersonic part of the pressure distribution dif- 
fers from the pressure distribution at the upper critical Mach number by 
nearly the same constant, although of opposite sign. 
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In order to investigate this difference further, the pressure dis- 
tributions at the upper and lower critical Mach numbers have been calcu- 
lated for the complete biconvex circular-arc airfoil and each of the four 
related airfoils having ms.y-innrm thickness forward and aft of the midchord 
station for which the results for Mach number 1 are shown in sketches (p) , 
(q), (r) , (s), and (t) . The results are shown in sketch (d*). It can be 



Sketch (d* ) 

seen that the three pressure distributions for each airfoil bear the same 
general relationship to each other as noted above, although the difference 
between the pressure distributions Is not always quite so constant as is 
observed for the wedge and circular-arc profiles. 


Pressure Drag 


Once the pressure distribution is known for a given airfoil, the 
pressure drag can be obtained directly by integration of equation (48) . 
The corresponding expression in terms of and Cp is 
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. _ [^y+l) 3 , 

'<3- T 5 / 3 ( 3 


2 


f 1 * a(z/t) 

J 0 ^a(x/c) 


a 



(84) 


Although the present theoretical results only permit the calculation 
of pressure drag for Mach numbers near 1 and for Mach numbers greater than 
the upper critical, these results, together with existing theoretical and 
experimental results, can be used to sketch the variation of pressure drag 
with Mach number throughout the transonic range. The airfoil for which 
the most information is available is, of course, the single-wedge profile 
for which an approximate solution for Mach numbers less than 1 has been 
given by Cole (ref. 36) and improved recently by Yoshihara (ref. 35 ) > 
that for Mach number 1 by Guderley and Yoshihara (ref. 32), and that for 
Mach numbers greater than 1 by Vincent! and Wagoner (ref. 34) , and experi- 
mental data have been given by Liepmann and Bryson (refs. 30 and 31 ) . 
Sketch (e’) shows a plot of these results, all in terms of the reduced 



parameters c^ and | OT and recast into the form consistent with the pre- 
sent formulation of the basic equations for transonic flow (see ref. 5, 

6, or 9 for additional information on this point), together with the 
results computed by use of the present theory. The new results are indi- 
cated by the solid lines, the former by dashed lines and by data points. 

The short vertical lines on the data points indicate Bryson’s estimate of 
the experimental accuracy of the data. As can be seen, the only point of 
difference between the present results and the previously existing results 
is at Mach numbers slightly in excess of the upper critical, and results 
from the error incurred in approximating the pressure jump through the 
bow shock wave by simple wave theory (i.e., by eq_. (15) rather than 
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eq. (6)). The positive slope of the drag curve at | M = 0 , or Mach number 1 
is in agreement with the result indicated by equation (51) and is typical 
for airfoils that do not close at the rear. 


Sketch (f 1 ) shows a summary of the comparable information for the 
front half of a biconvex circular-arc airfoil followed by a straight 
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section, for which experimental data have been given by Bryson in refer- 
ence 30. The theoretical values are again indicated by a solid line, 
and the experimental values by data points. 4 Although the amount of 
information available is much less than for the single-wedge airfoil, the 
results for both cases show striking similarity. 

Results for half airfoils are not typical of those for complete air- 
foils, however, as can be seen by comparison of the preceding results 
with the corresponding theoretical and experimental results illustrated 


4 The experimental values shown in sketch (f ! ) differ somewhat from 
those given originally in figure 21 of reference 30 because of the correc- 
tion of some inaccuracies in the calculation of the drag from the experi- 
mental pressure distributions given in figure 20 of reference 30. Although 
no explanation is known for the substantial negative drag indicated at 
subcritical Mach numbers and its existence must be indicative of some 
shortcomings of the experimental technique, its occurrence is an unmis- 
takable consequence of the measured pressure distribution. That this is 
so can be seen at a glance by comparison of the measured pressure distri- 
bution with that indicated by .linearized compressible flow theory, for 
which the drag is zero. 
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in sketch (g*) for conplete biconvex circular-arc airfoils. The experi- 
mental results are those given by Michel, Marchaud, and Le G all o in 



reference 37 and are obtained by integration of equation (48) together 
■with experimental values for the pressure distribution. The most promi- 
nent difference concerns the slope of the curve of c^ versus 
at ioo=0, or Mach number 1, for which the same procedures that led to 
positive values for a half airfoil, lead to zero slope for a complete 
airfoil. It can be seen that the experimental data support these values 
of the slope in both cases. Although the calculated values for drag are 
somewhat greater than those measured in the wind tunnel, most of the dis- 
crepancy can be attributed to the local effects of shock-wave boundary- 
layer interaction that occur near the trailing edge. Because this 
phenomenon depends on Reynolds number and may be of greatly diminished 
importance at full-scale conditions, Michel, Marchaud, and Le Gallo intro- 
duced, in the discussion of their experimental results, the concept of 
"extrapolated drag" to represent the drag that would occur in the absence 
of separation. This quantity is calculated by consideration of a pressure 
distribution that differs from the experimental pressure distribution in 
the vicinity of the trailing edge as a result of the replacement of the 
pressures actually measured by those obtained by extrapolation of the 
trends indicated at stations upstream of the separation point. Accord- 
ingly, the values for "extrapolated drag" given by Michel, Marchaud, and 
Le Gallo are also shown on sketch (g 1 ). As might be expected, the theo- 
retical values for drag are in better agreement with the values for 
"extrapolated drag" than with those obtained directly from the actual 
measurements . 

The results for biconvex circular-arc airfoils are typical of those 
for other complete airfoils. Attention is called, however, to the fact 
that the experimental values given by Michel, Marchaud, and Le Gallo in 
reference 39 for the airfoils that have the point of mavl mi im thickness 
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located forward of the midchord station are not bo reliable as those they 
give for airfoils that have a more rearward location of the point of max- 
imum thickness. This reduction in accuracy results from the facts that 
the method of testing and the fixed spacing of the orifices tend to dimin- 
ish the accuracy with which the contribution of the region near the lead- 
ing edge can be evaluated, and that the contribution of this region is, 
at the same time, of increased importance. 


Airfoils defined by 

Eq. (56) 

Eq. (61) 


Sketch (h* ) shows a 
summary of the calculated 
results for the pressure 
drag of the two families 
of airfoils described by 
equations (56) and (6l) 
with values for n rang- 
ing from 1.5 to 6. For 
this range of n the air- 
foils described by equa- 
tion (6l) have a range of 
location for the point of 
maximum thickness that 
extends from about 0.30 
to O.55 chord, and those 
described by equation (56) 
have a range of locations 
extending from about 0 . 4-5 
to O.7O chord. In addi- 
tion to lines for con- 
stant which correspond 

to lines of constant Mach 
number for a group of air- 
foils having the same 
thickness ratio, a line is 
also shown for loo cr > which 
corresponds to the line for 
the lowest Mach number for 
which the flow is purely 
supersonic. It can be seen 
that the variation of pres- 
sure drag with the location of the point of max-l mum thickness at Mach 
number 1 is quite different from that indicated by lines of constant £ 
or Mach number, for purely supersonic flows, but is rather similar to 
that indicated by the line for g . An interesting feature of the 
results for purely supersonic flow is that the drag is not the same in 
forward and reverse flow, as is indicated by linearized compressible 
flow theory (see ref. 14-3 44 ). 
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CONCLUDING REMARKS 


It appears worthwhile, in conclusion, to summarize and contrast the 
alternative discussions presented in the main text and in the Appendix of 
the general procedures involved in the approximate solution of all the 
problems treated in this paper. The arguments presented in the Appendix 
are based essentially on the idea of diminishing the importance of the 
higher order terms, and hence concentrate on the contributions stemming 
from the double integral of the integral equation appropriate for each 
case. The arguments presented in the main text lead to the same conclu- 
sions, but are based essenti ally on the idea of linearizing the transonic 
equation in a small region by replacing part of the nonlinear term by a 
constant A, and then introducing different values for A for different 
points in the field. This procedure might be considered equivalent, in 
some sense, to the replacement of the original nonlinear equation by a 
different linear differential equation for each point. Results obtained 
by solution of the equations at this stage depend, of course, on the 
choice of A and must be assembled in order to determine the final 
results. This step is accomplished in each case by putting the results 
into such a form that a first-order nonlinear ordinary differential equa- 
tion is obtained upon substituting for A the quantity it originally 
replaced. At this point, the equations encountered coincide with those 
obtained following the procedures described in the Appendix and the 
remainder of the analysis proceeds in identical manner. In the cases con- 
sidered herein, the differential equation is always of sufficiently simple 
form that it can be integrated analytically and the result expressed in 
closed analytic form. This integration Implicitly introduces the assump- 
tion of continuity of the velocity or pressure distribution and leads, 
upon evaluation of a single constant of integration, to the final result. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Mar. 11, 1957 
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APPENDIX A 


REINSPECTION AND ANALYSIS OF METHOD OF APPROXIMATION 


INTRODUCTION 


The methods used In the main text of this paper are appealing for 
their brevity and for the efficiency with which approximate solutions of 
the nonlinear* equations of transonic flow theory are found. These methods 
are not entirely satisfying to the critical reader, however, because cer- 
tain elements appear to be arbitrary and there is no a priori way in which 
the accuracy of the approximations can be judged. Both the virtue and 
weakness of these procedures are the result of introducing the essential 
simplifications at the beginning of the analysis. If the introduction of 
approximations is deferred to the end, the relations that occur in the 
initial stages of the analysis are, of course, more complicated than those 
presented in the main text. Consideration of these relations is, however, 
essential for an understanding of the basis for the method of 
approximation . 

The following discussion of the approximate solution of problems of 
transonic flow theory is based on consideration of integral equations 
derived from the differential equations of transonic flow theory by stand- 
ard application of Green 1 s theorem. Since the details of each of the 
three cases, hyperbolic, elliptic, and parabolic, are somewhat different, 
each case is considered separately. In each case, exact relations are 
retained as long as possible and the approximations, when introduced, are 
seen to be closely related to those employed in the method of successive 
approximations commonly used in the determination of higher approximations 
in compressible flow theory. The following paragraphs will be concerned 
at the outset with the derivation of integral equations for each of the 
three cases, and subsequently with the discussion of the simplification 
and approximate solution of these equations. 


DERIVATION OF INTEGRAL EQUATIONS 


All of the subsequent analysis proceeds from Green’s theorem. There 
are many forms of Green 1 s theorem, but a sufficiently general form for 
all of the present purposes is that associated with the linear opera- 
tor L(g) defined as follows (ref. 45 f pp. 244 - 247 ) : 


L( G) = AD^qc 


(Al) 
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where A and B are constants. Green’s theorem states that the follow- 
ing relation holds between any two arbitrary functions Q and having 
continuous first and second derivatives : 


JJ +L(fl)-fiM(+) an 



z 



in which R refers to the interior_of an arbi- 
trary region bounded by the curve C, as shown 
in sketch (i’), M is called the adjoint 
differential operator 

M(t) = Bxx + * zz - Bf x 
A represents the quantity 



and b/bv stands for a derivative in the direc- 
tion v, and can be written as a linear 
differential operator. 

^ - f2xCOs(v,x) + G z cos(v,z) 


The direction v is called the conormal, and its direction cosines are 
related to those of the normal n according to 


A cos(v,x) = A cos(n,x) 
A cos(v,z) = cos(n,z) 


Hyperbolic Case 


The initial step in the present derivation of the integral equation 
appropriate for the discussion of the hyperbolic case is to subtract 

from both sides of equation ( 3 ) , and to write the resulting equa- 
tion as follows: 


-^H9xx + <P ZZ = (^-l+ku-^)^ = f H (A3) 

The symbol Ag refers to any positive finite constant. The form of 
Green’s theorem associated with the linear operator 


(Afc) 


L(q>) = Lh(cp) = -Ah'Pxx + ^zz 
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will now be applied , whence 

A = -Ag , 3 = 0, %(cp) = L H (cp) (A5) 

The quantity £2 is now identified with the perturbation potential cp, 
and i|r with an elementary solution of Mjj(i|r)=0, in particular with the 
unit supersonic source cr^ defined as follows: 

for (x-i) 2 >A H (z-p 2 (A6a) 

2-JAh 

0 for (x-|) 2 <Ag(z-£) 2 (A6b) 

If the region 5 is so selected that the inequality expressed in equa- 
tion (A6a) holds throughout, equation (A2) reduces to 


- J % It 113 = JJ = JJ f H® ( A 7) 

® H R R 

in which Ag- and v^- refer to the special forms of A and v consistent 
with equation (Ak) or (A5) , the running coordinates of integration 

|,£, and the field point at which (p is to be evaluated is x,z. 
Equation (A7) is now applied to the region indicated in sketch (j*). 
Note that the wing, wake, and shock waves must 
be excluded from the region of integration. 

It should be noted that sketch ( j 1 ) is only a 
schematic illustration to help define the 
quantities involved in equation (A 7 ) , and that 
the shock wave, indicated as a detached bow 
wave, might instead be attached to the leading 
edge, to the trailing edge, or situated some- 
where along the chord. The single shock wave 
illustrated in sketch (j*) could, moreover, be 
replaced by a campli£ated system of shock 
waves. The region Rjj may likewise be con- 
sidered to consist of a single region, or as 
the sum of a number of regions as typified by 
adding to Rjj the region situated upstream of 
the bow wave and within the dotted lines of Sketch (j*) 

sketch (j*). In any case, the general considerations remain the same, 
and the following integral equation is obtained for cp at an arbitrary 
point x, z : 
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<p(x,z) = 


1 

2^ 







1 

2 n/^B 





1 

2-/^H 


Rs 


(A8) 


where 


d<P _ 52a 52i 

^ = X ~ ZZ 


Acp = cp u - cp z 


in which u and l refer to conditions on the upper and lower sides of 
the wing and wake, where 

8 59 5cp a ^ 5cp b 
Sv H ” 5v H + Sv H 

in which a and t refer to conditions on the two sides of the shock 
wave Sg, and where 


a<p _ 

a v 


■>H cos(n,x) 


acp / ^ 

+ — cos(n, z) 

35 


Elliptic Case 


The integral equation appropriate for the discussion of the elliptic 
case can he derived by use of procedures analogous to those described in 
the preceding paragraph for the hyperbolic case. The initial step is to 
add AgCpxx to both sides of equation (3), and to write the resulting 
equation as follows: 


VPjcx + ^zz 


A E -(l-*k 2 -ku) 

- 


?xx = f E 


(A 9) 


As in the hyperbolic case, the symbol Ag refers to any positive finite 
constant. The form of Green* s theorem associated with the linear operator 
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(A10) 


(All) 

The quantity ft is now identified with cp, and i(r with an elementary 
solution of Mg(i]f)= 0 , in particular the unit subsonic source defined by 
the function 

t - ® E i= in J (x-D^Cz-O 2 ("2) 

In this way equation (A2) reduces to 


L(cp) = Lg(cp) = AgC p^ + $ zz 

will now be applied, whence 

A = Ag , B = 0 , Mg(cp) = LgCcp) 



“ JJ 


R 


CT E f E dR 


(A13) 


in which Ag and v E refer to the spe- 
cial forms of A and v consistent with 
equation (AlO) or (All) , the running 
coordinates of integration are 
and the field point at which 9 is to 
be determined is x,z. If equa- 
tion (A13) is now applied to the 
region r E surrounding the wing, wake, 
and shock waves, as Illustrated in 
sketch (k* ) , and the a priori assump- 
tion Is made that the perturbation 
field attenuates sufficiently fast with 
distance to negate the contribution of 
the surface integral over the large 
circle in the limit as the radius goes 
to infinity, the following integral 
equation is obtained for cp at an 
arbitrary point x,z: 



Sketch (k 1 ) 




— Am — 

°RA;e 5 - — dSg + 

Sv E 


JJ 


crgfgdRg 


*E 


(AU) 


where A 9 and A(d9/d£ ) have the same meaning as in equation (AS), but 
where 
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B d? d<P a } 

^E ^E ^E 

in which a and b refer to conditions on the two sides of the shock 
wave Sp, and where 


Ae j-L = A e cos(n,x) 
OVg 85 


+ — cos(n,z) 


Parabolic Case 


The initial step in the derivation of the integral equations appro- 
priate for the discussion of the parabolic case is to subtract Apcp x 
from both sides of equation (3) , and to write the resulting equation as 
follows : 

^zz-VPx = (Mnc 2 - 1 )T xx +(kcp xx -?v p )cp x = fp (A15) 


The symbol Ap refers again to a finite constant, which may be either 
positive or negative. The form of Green's theorem associated with the 
linear operator 


L(<p) = Lp(cp) = 9 zz -ApCp x (A16) 

will now be applied, whence 

A = 0 , B = -A P , Mp(cp) = cp zz + A p (p x (A 17 ) 

The quantity is again identified with 9, and t with an elementary 
solution of Mp(\|f)= 0 , in particular with the function 


r 



Ap(z-e) 2 ' 

4 (x-|) . 


0 


Ap ^ 

for > 0 

x-s 


Ap 

for < 0 

x-6 


(Al8a) 

(Al8b) 




This function assumes a role in the analysis of the parabolic case that 
is analogous to that of the unit subsonic and supersonic sources in the 
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e lli ptic and hyperbolic cases. _In mathematical literature, the linear 
partial differential equation Lp(cp)=0 with positive Ap arises in the 
study of heat conduction, and the function cfp is often referred to as 
a unit heat source. If the region R is so selected that the inequality 
expressed in equation (Al8a) holds throughout, equation (A2) reduces to 


- J Ap (a-p -q> ^ “TvpCTptp cos(n,x) ds = JJ OpIp(cp)dR = JJ' CpfpdR 


R 


(A19) 


in which Ap and vp refer to the special forms of A and v consistent 
with equation (Al6) or (A17), the running coordinates of integration 
are £,£, and the field point at which is to be determined is x,z. 

It is apparent from the condition imposed on Rp, that Ap/(x-|) is 
greater than or equal to zero, that two distinct subcases result depending 
on the sign of Ap. 


Positive Ap.- If Ap is positive, the 

region Rp appropriate for the application of 
equation (A19) is that part of space upstream 
from the point x,z. Again the wing, wake, 
and shock must be excluded from the region of 
integration, as illustrated schematically in 
sketch (z 1 ). If equation (A19) is now applied 
to the region Rp, and it is assumed that the 
perturbation field attenuates sufficiently 
fast with distance to negate the contribution 
of the surface integral over the outer bound- 
ary in the limit as the radius is increased 
to infinity, the following integral equation 
is obtained for cp at an arbitrary point x, z: 


/ 

/ Rp 


UcD 


•x,z 




gp-' 


Cp'i 


Sketch ( Z * ) 


cp(x,z) = 




dpA — - A<p 


W 


) il "k£ a ^ (I7) 

Sp Rp 

(A20) 


where A9 and A(d<p/d£) have the same meaning as in equation (A8) , but 
where 


= ^a + 

Svp Svp Svp 
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in which a and b refer to conditions on the two sides of the shock 
wave Sp, and where 


-Sep Sep 
* = a; °° s(n ' z> 


U (0 



Sketch (m 1 ) 


Negative Ap.- If Ap is negative, 

the region Rp appropriate for the appli- 
cation of equation (A19) is that part of 
space downstream from the point x,z. If 
equation (A19) is applied to the region Rp 
surrounding the appropriate part of the 
wing, wake, and shock waves, as illustrated 
in sketch (m* ) , and the contribution of the 
surface integral over the outer boundary 
vanishes as the radius is increased to 
infinity, the following integral equation 
is obtained for cp at an arbitrary 
point x, z : 


cp(x,z) 



+ -r- rr 


Ap JJ 

R P 


(A21) 


where the symbols have the same meaning as in equation (A20) , except 
that Sp and Rp now refer to those portions of the shock waves and space 
situated downstream from x,z. 


SOME PROPERTIES OF THE INTEGRAL EQUATIONS 


Although the four integral equations derived in the preceding para- 
graphs and written explicitly in equations (A8) , (Alt) , (A20) , and (A21) 
axe quite different in most respects, they do possess a number of prop- 
erties in common that are of concern in the present discussion. Perhaps 
the most obvious similarity is that each integral equation consists of a 
term that involves integration over the wing and wake, another term that 
involves integration over the shock waves, and a third term that involves 
integration over the surrounding space. The integrals extend over all 
space in the elliptic case, but only over part of space in the parabolic 
and hyperbolic cases. It is important to realize that there is no direct 
connection between the region of integration in each of the integral equa- 
tions and the region of dependence in the solution, or in the physical 
flow, and that these two regions may, in fact, be distinctly different 
in some applications. 
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The first term, in each, integral equation involves a distribution of 
sources a proportional to A(S<p/&0 and doublets bc/bt, proportional 
to A?. Since A(dcp/c)£) is equal, according to the boundary condition 
given in equation (l) , to U^ACcLZ/dx) and AP is proportional to an 
x-wise integration of the lift, it follows that the part of the term con- 
taining sources is associated with the thickness distribution, and that 
part containing doublets, with the effects of camber and angle of attack.. 
The latter part of the first term is zero for al 3 of the nonlifting air- 
foils discussed in the main text of this paper. The first term in the 
elliptic and hyperbolic cases is familiar in compressible flow theory 
because, upon equating Xg to l-l^ 2 or Xg to l^-l, it reduces to the 
well-known solution for cp in the linearized theory of subsonic and 
supersonic flow around thin airfoils. 


The second term in each integral equation involves a distribution of 
sources a proportional to AS(d<p/dv) over any part of the shock waves 
that is situated in what otherwise would be part of the region of inte- 
gration R. There is no doublet distribution on the shock wave, as on 
the wing and wake, because cp is continuous across the shock wave. 
Although the contribution of the integral over the shock waves is often 
difficult to evaluate because 8(S<P/Sv) is unknown and must be determined 
as part of the solution, there are a number of important applications in 
which this term either vanishes completely, or contributes nothing to the 
values for cp along the chord of the airfoil. The simplest class of 
problems for which this term vanishes is, of course, that in which the 
flow is subsonic everywhere and is hence shock- free. The contribution of 
this term will also vanish in parts of the field even if shock waves are 
present, provided they are situated entirely downstream of the region of 
integration in the hyperbolic case or the parabolic case with posi- 
tive Xp, or entirely upstream of the region of integration in the para- 
bolic case with negative Xp. The contributions of the shock wave vanish 
in the above situations because the complete term disappears from the 
integral equation. If the term remains, however, each element of the 
shock wave provides a contribution to cp that depends upon its str eng th 
and orientation. There are, moreover, certain directions in which an 
element of a shock wave can be oriented that result in no contribution 
to cp in the parabolic and hyperbolic cases. Thus, in the parabolic 
case, the contribution vanishes when the element of the shock wave extends 
perpendicular to the x axis, so that cos(n,z) in the second term of 
equation (A20) or (A21) Is zero. It Is s imi larly evident from equa- 
tion (A8) that an element of a shock wave contributes not hing to cp in 
the hyperbolic case if the direction cosines and velocity components on 
the two sides of the shock wave satisfy the relationship 



cos(n b ,x) 


(A22) 
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This relation can be expressed completely in geometric terms by combining 
it with the following equation that can be derived from consideration of 
the fact that cp is continuous across the shock 

cos (n-jj j z) ^ w a -Wfr 
cos(n v x) Ua-u-b 

Thus, the contribution to cp vanishes if the element of the shock wave 
is oriented bo that 


COBpbjZ) 

cos(n-b,x) 


2 


tan 2 (n- b ,x) = A s 


(A24) 


If the same lin e of reasoning is applied to the elliptic case, relations 
analogous to equations (A22) through (A24) occur in which Ah is replaced 
by -Ag. From such considerations it would appear that the contribution 
to cp vanishes if the element of the shock wave is oriented so that 


COB ( n b* z ) 

cos(nb,x) 


2 


= 


(A25) 


Since A E is required to be a positive quantity, however, it is clear 
that there is no orientation for which the contribution vanishes. In the 
hyperbolic case, on the other hand, A s is required to be positive and 
there are always two particular orientations for which an element of a 
shock wave contributes nothing to 9 . It is interesting to note, before 
leaving the discussion of the second term of each of the integral equa- 
tions, that the particular orientation for which an element of a shock 
wave provides no contribution to cp coincides in all three cases, ellip- 
tic, parabolic, and hyperbolic, with the directions of the characteristic 
lines of the associated form of the linear partial differential equation, 
L(cp)=0. The reader should observe, however, that these characteristic 
lines have no particular physical significance, inasmuch as their exist- 
ence and direction depend on an arbitrary choice of a value for A. 

The third term in each integral equation involves integration of the 
effects of a distribution of sources a proportional to f over that 
part of space surrounding the airfoil that is enclosed within the region 
of integration R. The contribution of this term does not vanish, except 
in almost trivial circumstances such as occur, for example, in the hyper- 
bolic case for points upstream of a bow wave provided A s is equated 
to Moo 2 -! so that % is zero. Discussion of the contribution of the 
third term in each integral equation will consequently constitute the 
subject of much of the remainder of this Appendix. An interesting 
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property of each of the Integral equations that is worth noting before 
proceeding to the more specific discussion of each case is that the inte- 
grated strength of the sources in all space exterior to the airfoil, 
including those distributed along the shock waves, is equal, but opposite 
in sign, to the integrated strength of the sources distributed along the 
entire chord of the airfoil. 


SALIFICATION AND APPROXIMATE SOLUTION OP 
THE INTEGRAL EQUATIONS 


No general methods are known for the analytical solution of the 
integral equations given in equations (A8) , (Alh) , (A20) , and (A21) . 
Although certain simplifications can be made by restricting attention to 
nonlifting cases, and to cases in which the shock waves are not in the 
region of dependence, the essential difficulties remain because the inte- 
gral equations are nonlinear, just as is the differential equation from 
which they are derived. The principal method that has been employed in 
the past for the solution of similar problems is that of successive 
appro xi mation in which cp is expanded in a power series of some param- 
eter such as the thickness ratio, and the coefficients in this serieB are 
determined as the solution develops. The first approximation in these 
methods is generally either the solution for linearized compressible flow 
theory or for incompressible flow, and the second and higher approxima- 
tions are determined by iteration procedures in which linear equations 
are solved at every step. In practice, these methods have been found very 
difficult to apply to problems of compressible flow, and calculations of 
higher approximations than the second have, in most cases, proved prohib- 
itively lengthy. (See ref. 3 for a resume.) Serious questions of con- 
vergence remain in the existing solutions of this type, and it is doubtful 
if the results apply when mixed subsonic- supersonic flow occurs. 

Another type of successive approximation procedure which involves 
the solution of quadratic equations with every iteration step is des- 
cribed in reference 9 and applied to the calculation of pressure distri- 
butions on circular-arc airfoils for al 1 Mach numbers up to unity. 

Although the calculations could only be accomplished after the introduc- 
tion of certain approximations , whose influence on the result is difficult 
to ascertain, the general procedure appears to succeed even with the 
occurrence of mixed flow. 

The methods applied in the main text of this paper can be considered 
as the first step of still another type of successive approximation pro- 
cedure in which certain nonlinear features are incorporated into even the 
first approximation . This procedure possesses the advantage of yielding 
results that disclose much of the nonlinear effects in the first approxi- 
mation, and of making unnecessary, in many applications, the difficult 
task of iteration to determine higher approximations. A simple heuristic 
description of the analysis is given in the main text. The following 
paragraphs are concerned with a more detailed examination of the 
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approximations involved in the analysis, and of the relationship between 
the present approximation and the approximations previously employed in 
the solution of problems of compressible flow. Although it is apparent 
that much of the discussion could be applied to lifting airfoilB, atten- 
tion is confined, as in the main text, to symmetrical nonlifting airfoils 
for which the following relations are to be applied in the first term of 
each integral equation: 


A9 = 0 , 



2U„ 


dZ 

d| 


Hyperbolic Case 


(A26) 


The first problem to be discussed is the approximate solution of the 
integral equation given in equation (A 8) for the hyperbolic case under 
the restrictions that the flow is purely supersonic so that 

l^o 2 -l+ku>0 (A27) 

Equation (A8) contains both a line integral over the shock wave Sg and 
a double integral over the surrounding region Rg. The double integral 
can be integrated x-wise, however, because the integrand is a perfect 
differential. This partial integration results in a term that exactly 
cancels the line integral on the shock wave and equation (A8) can be 
rewritten as follows: 


cp(x,z) 



The symbol > below the integral sign of the second t erm indicates that 
the integration is to be carried over the lines j*s=x- Ag( | z-£ |) extending 
upstream from x,z. Differentiation yields the following relations 
for u and bu/bx: 


u(x,z) 


(A29) 

-Ik- Z" (x-JtJ \z\) - — -f U (moo 2 -1+§: -Ah) dC 

(A30) 


du(x,z) 

bx 
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It should be noted that equations (A28) , (A29) , and (A30) are all integral 
equations and that each is an exact relation valid for any positive value 
of Ag. 

The results of linearized supersonic flow theory, as well as those 
of the counterpart, in the sma ll disturbance theory of transonic flow, of 
higher approximations can be reproduced from any of the above equations 
upon equating Ah to 1^ 2 -1. In this way, all contributions of first order 
are included in the first term, on the right and the contribution of the 
second term on the right is, at most, of second order. Thus, the familiar 
expression of linearized supersonic flow theory follows immediately upon 
disregarding the contribution of the second term 

u(x,z) U " Z' + °(u 2 ) (A31) 

Its counterpart, correct to second order in Z* , can be determined by 
application of the method of successive approximations in which the con- 
tribution of the second term is approximated by replacing the unknown u 
by the first-order approximation provided by equation (A3l) and integrat- 
ing. In this way, the following result is obtained for points on the 
airfoil. 


u(x,0) 


Uqc dZ _ kUco g ( dZ Y+ 0 ( 62 
** " 4(I^ 2 -1) 2 V A* / V dx 


(A32) 


Note that although the first- and second-order approximations for u are 
different in general, first-order theory is sufficient to determine to 
second-order accuracy the point where u vanishes, that is where dZ/dx=0. 
First-order theory, moreover, provides the exact location for the point of 
zero u in the absence of shock waves, but this simple result is true 
only to second order if there are shock waves situated within the Mach 
forecone of the point (x,0) . Although the difficulties of integration 
are such that only the first few steps of the method of successive approx- 
imation can be evaluated in all but the simplest examples, the method can, 
in principle at least, be repeated indefinitely to establish the result 
accurate to any desired order. The result for u on the airfoil surface 
appears in the form of a power series involving ascending powers of 

(dZ/dx)/(jf^ 0 2 -l) 3 ^ 2 . It is clear that a valid approximation is obtained 
provided the absolute value of dZ/dx is sufficiently sma] J at all 
points, and the Mach number is not too close to unity. The failure 
associated with excessive positive values for dZ/dx usually occurs near 
the leading edge of the airfoil and is associated physically with detach- 
ment of the bow wave and the occurrence of local regions of subsonic flow. 
The failure associated with excessive negative values for dZ/dx, and 
clearly illustrated in sketch (c) of the main text, usually occurs near 
the trailing edge of the airfoil and is purely mathematical in origin. 
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In actual practice where only the first term or two may he evaluated, the 
result fails to provide adequate information regarding the ultimate con- 
vergence or divergence of the series and the question must he settled in 
each application hy appeal to more exact solutions. It is important to 
realize that these uncertainties are not inherent in the integral equa- 
tions given in equations (A28) through (A30) , hut enter the analysis with 
the assumption that the solution can he approximated satisfactorily hy 
application of the particular form of the method of successive 
approximation described above. 

The method of analysis employed throughout the main text iB equiva- 
lent, from the present point of view, to the first step of a slightly 
differ en t method of successive approx im ation that proceeds from consider- 
ation of the infinite set of relations that result if different values are 
selected for Ah in the determination of conditions at different points 
in space. Analyses based on such a system of equations are more complex, 
in general, than those based on a single equation, hut this increase in 
complexity is counterbalanced, in the present applications, hy the fact 
that approximations can he In troduced on the basis of local, rather than 
global, considerations. This fact makes possible the incorporation of 
some of the higher order or nonlinear contributions, as well as all of 
the first-order contributions, into the first term on the right in each 
relation of the infinite set, thereby reducing the contribution of the 
term containing the unknown u( g , £ ) . If a rule for the selection of Ah 
can be found that achieves this effect and if it can be expressed in ana- 
lytic form, the infinite set of relations can be expressed once again In 
the form of a single equation; and the remainder of the analysis can pro- 
ceed in a manner analogous to that described in the preceding paragraph 
for the classical method of successive approximation. 

The method employed in the analysis of the hyperbolic case in the 
main text of this paper is equivalent to the first step of a successive 
approximation procedure based on the infinite set of equations for 5u/6x 
typified by equation (A30) with Ah equated to the local value of 

-l+ku(x,z) in each relation. Although it would appear from cursory 
examination of equation (A30) that the error incurred in so doing would 
be of second order, it will be demonstrated below that the result is 
actually accurate to second order and that the error is, at most, of third 
order. Although equations (A29) and (A30) , upon which the remainder of 
the discussion, of the hyperbolic case Is based, are exact within the 
approximations of transonic flow theory, they are not in the most advan- 
tageous form for the following discussion because of the presence of the 
integral in the second term of each relation. If, however, attention is 
confined to the evaluation of the result at the airfoil surface to second- 
order accuracy, and to cases for which Ah is restricted to values that 
differ from Moo-l "by* at most, a quantity proportional to u(x,z), that 
Is AH=H» £- i+ au ( x ^ z ) > the integrals can be integrated and the following 
relations result: 
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u(x, 0 ) 


Uco dZ 

■T*** 


u -*h) 

™ 2 (^- 1 ) 


+ 0 (u 3 ) 


(A33) 


5u(x, 0) _ _ Uco d g z _ ^xCHx^-i+ku-Ag) + q(u s ) 

^ ^ 2 2 ( 1 ^ 2 - 1 ) 1 " U ' 


(A3^) 


It is now clear that the substitution of M^-l+ku for Ag in equa- 
tion (A34) yields 


du(x, 0 ) 

dx 


— — S + 0 (u=) 
, 2 -l+ku ^ 


(A35) 


This result corresponds to equation (ll) of the main text and leads, upon 
integration and insertion of the auxiliary relation that u vanishes 
where dZ/dx is zero, to equation ( 15 ) relating Cp and dZ/dx. This 
relationship is commonly designated as that of simple wave theory. From 
the above discussion, it is clear that equation ( 15 ) must be correct to 
at least second order, as indeed simple wave theory is known to be for 
the pressures on the surface of an airfoil. If there are no shock waves 
in the region influencing conditions at the point (x, 0 ), the flow field 
is characterized by a single family of waves; and it can be shown that 
the error term, indicated in equation (A35) vanishes completely. The 
resulting relation Is thus exact within the approximation of transonic 
flow theory. It is interesting to observe that the use of the same 
relation for Ag in equation (A29) for u leads to 


H 22 + ku 2 + 0(uS) 

A, 2 - ^ 


(A36) 


and results in errors of second order if only the first term is used. 

There is another choice for Ag that is not mentioned in the main 
text that will remove the second-order error if only the first term of 
equation (A29) or (A33) for u is used and that is Ag=^ ;o 2 -l+(ku/ 2 ) , 
since then 
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u 


U Q 



dZ 

dx 


+ 0(u 3 ) 


(A37) 


This result is not the same as that of simple wave theory, hut is 
recognized upon rearranging so that 

u iTii = -u*, g (A 38 ) 

as being the square root, with appropriate choice of sign, of the shock 
relation given in equation (6) with Ug, and w a equated to zero. The 
result obtained by use of equation (A 38 ) is thus equivalent to that 
obtained by equating the pressure at each point of the airfoil surface 
to that on a tangent wedge. Such a procedure has been proposed previously 
and is sometimes called ” tangent -wedge theory." Although the first two 
terms of the formal expansion of either equation (A 38 ) of tangent-wedge 
theory or equation (l4) of simple wave theory agree with the second-order 
result obtained by use of the method of successive approximations and 
given in equation (A 32 ) , the results of either simple wave theory or 
tangent -wedge theory are to be preferred in applications because they 
approximate the proper termination of the solution when dZ/dx becomes 
too large, and do not fail spuriously at larger negative values of dZ/dx. 


Elliptic Case 


The second problem to be discussed is the approximate solution of 
the integral equation given in equation (Al4) under the restriction that 
the flow is purely subsonic so that 

l-Moo 2 -ku> 0 (A39) 

This restriction implies that the integral over Se is zero, since there 
are no shock waves in a purely subsonic flow. Thus, equation (Al4) reduces 
to 

c +0 ° 

<p(x,z) = 2 U m J ff E if d| + JJ f E cr E d| d£ (A40) 


Differentiation yields the following relations for u and 8u/8x 
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u(x,z) 



Scr E dZ 
— Si — d| 
Sx dg 



-OO 


(Aia) 


chi(x,z) 

dx 


2Uco 



22 ai + 

dx 2 dg 


IT 

-00 


f- 


E 



d| dg 


(A42) 


where 

ff E ~= in N /(x-|) 2 + A E (z-0 2 

2jr^ 


^ a E 

Sx 



*-s 

(x-g) 2 +A E (z-g) 2 



Bx 2 



(x-g) 2 -A E (z-g) 2 

^ 2 I 2 

(x-g) +A E (z-g) 

- _ 


f E “ 


A E -[l-l^ 2 -ku(g,g) 


Su 


Equations (AtO) , (A4l) , and (Alf-2) are all integral equations and each is 
exact for any positive value for A E . The solution of any of these inte- 
gral equations is complicated not only because the relations are non- 
linear, but also because the kernel, designated by cr E or its derivative, 
is infinite at the point g=x, g = z. 


The familiar result of linearized subsonic flow theory can be 
obtained from any of the above equations by replacing A E by 1-M^ 2 a nd 
disregarding the contribution of the double integral as being of higher 
order. The corresponding result, correct to second order, can be 
obtained by application of the method of successive approximations in 
which the contribution of the double integral is approximated by means 
of the first-order result to evaluate f E at each point. Higher approx- 
imations can be obtained, at least in principle, by repeated application 
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of the same procedure, except that fg is evaluated at each step by use 
of the results of the next lower approximation. In this way, an approxi- 
mate expression for. the solution is determined in the form of a truncated 
power series. Although the difficulties of integration are so great that 
few cases have been evaluated beyond the second approximation, it appears 
that the process converges to the desired solution for thin airfoils pro- 
vided that certain well-known difficulties associated with stagnation 
points are properly accounted for and that, as in the hyperbolic case, 
the Mach number is not too close to unity. Again the results provided by 
the method of successive approximation indicate no definite l imi t for the 
Mach number. Comparison with experimental results shows that the trends 
displayed by the results are generally confirmed for Mach numbers less 
than the critical Mach number, but are essentially refuted for greater 
Mach numbers. 

It is interesting and informative to compare the results obtained in 
the manner described above with those obtained by application of an alter- 
native version of the method of successive approximations described in 
reference 9 that involves the solution of quadratic rather than linear 
algebraic equations at each step of the iteration process. The equation 
fundamental to this discussion Is obtained from equation (A4l) by again 
equating to l-M^ 2 and integrating the double integral by parts. In 
this way the following integral equation is determined for u: 


u - 


u 2 




= 2U(v 


+00 

as-* ff 

J 0 dx d£ JJ 


U 2 (g,0 & * E 

2 5x 2 


*6 (A43) 


Although equation (A4-3) is completely equivalent to equation (A4l) , it 
is, in certain respects, superior from the point of view of obtaining 
approximate solutions. This is because the predominant effects of the 
region near the point x,z, which form a major contribution to the value 
of the integral in equation (A4l) , are furnished in equation (A43) , by 
the term involving the square of u standing outside the integral. 
Although the difficulties of integration are as great or greater than 
encountered in the classical method of successive approximations and only 
the first few steps can be evaluated without approximation in any specific 
application, certain general features of the solution are clearly defined. 
In particular, it is shown in the report version of reference 9 that the 
results obtained for flows that are subsonic everywhere converge, in the 
limit of an infinite number of iteration steps, to the same result as 
ultimately obtained by application of the classical method of successive 
approximation. Whereas there is considerable doubt about the precise 
range of convergence of the iatter result, the result obtained by appli- 
cation of the quadratic method of successive approximation clearly termi- 
nates with the occurrence of sonic velocity somewhere in the flow. The 
termination of the solution is recognized by the disappearance of real 
roots of the quadratic equation and is apparent at every step of the 
iteration process. It is evident, moreover, from comparison of the two 
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sets of results that the series expansion for the solution obtained by 
the classical method of successive approximation converges only for purely 
subsonic flows and that the results indicated for mixed or transonic flows 
are false. These properties of the result obtained by the quadratic 
method of successive approximation are consistent with the numerous argu- 
ments and proofs for the nearly nonexistence of continuous shock-free 
transonic flows that have been advanced in recent years. (See ref. 46 
for a brief resume.) 

This difference in behavior can be readily Illustrated if one con- 
siders the expressions for the pressure coefficient at the midpoint of a 
symmetrical circular-arc airfoil that are obtained following the comple- 
tion of the first two steps of each iteration process. The result pro- 
vided by use of the classical method can be readily obtained from the 
third-order result quoted in equation (37) of the main text and is 


Cp 


8 

jt 


sIl-Mo? 


{«-$) 


( y +1)]^ 2 


T 2 + 0(t s ) 


(A44) 


The first term represents the result obtained if one considers equa- 
tion (A4l) and disregards completely the contribution of the double inte- 
gral. The result is precisely that of linearized compressible flow theory 
for the particular point under discussion. The second term represents an 
approximation for the contribution of the double integral obtained by 
replacing fg for each point in space with the result provided by line- 
arized theory. The results obtained by application of equation (A44) 
Indicate that -Cp increases indefinitely with increasing value of 
t/(1-H» 2 ) 3/2 and appear to apply for mixed, or transonic, flows as well 
as for purely subsonic flows. The corresponding results obtained by use 
of the quadratic method of successive approx? mat ion are found by consid- 
eration of equation (A43) • A first approximation obtained by disregarding 
completely the contribution of the double integral is 


2(l-^ g ) 

*to 2 (7+l) 


L 8 ^ a (r~ 
J " * (i-iO s/2 


t+0(t 2 ) 


(A45) 


and a second approximation obtained from the evaluation of the contribu- 
tion of the double integral by use of the result provided by lin earized 
theory for f-g at each point is 
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2(l-^ g ) 

*k 2 (r+i) 




8 ^(7+1) T+ (± + i\ *k 4 (r+i ) 2 
51 d -^ £ ) 3/2 V * 2 2 ) ( 1-O 3 


t 2 +0(t 3 ) 


(A46) 


The results obtained by application of either equation (A45) or (A46) also 
indicate that -Cp increases with increasing values of t/(1-Mso 2 ) 3/ ^ 2 , 
but that this trend terminates when Cp reaches the critical value asso- 
ciated with the occurrence of sonic velocity, that is, when 


_ _ 2 ( 1-0 _ « 

05 " Or+i) “ ° Pcr 


(5) 


At this point, the slope of a curve representing the variation of G p 
with Mqo Is infinite. In spite of these distinct differences in behavior, 
it is important to note that a formal power series expansion of equa- 
tion (A46) in terms of t agrees to second order with equation (A44) , 
and that this agreement increases by one order of t upon the completion 
of each additional iteration step. Results obtained following completion 
of additional Iteration steps continue to follow the same trends. Those 
indicated by the classical method never provide any information regarding 
the precise range of convergence, and those indicated by the quadratic 
method always terminate with the occurrence of sonic velocity. 

The method employed In the main text can be considered as an alter- 
native procedure devised in an attempt to improve the quality of the first 
approximation and to diminish thereby some of the necessity for the eval- 
uation of higher approximations . Before proceeding, it is important to 
recall that equations (A40), (A4l) , and (A42) are all integral equations 
valid for any positive value for Ag; and. that each can be considered, In 
the same way as described for the hyperbolic case, as a typical member of 
an infinite set of relations that result if different values are selected 
for Ae In the determination of conditions at different points in space. 

It appears plausible that an increase in the accuracy, although not the 
mathematical order, of the first approximation might occur If Ag Is 
equated not to simply 1-M^ 2 , but to l-M aj s -ku(x,z) because then the 
function fg In the double integral of each integral equation reduces 
to k[u(x,z)-u(|,(;) ] (8u/d|) and hence vanishes at the point |=x, £«=z 
where cr E , ciffg/Bx, or 8 2 o '-g/dx 2 are infinite. If this procedure is 
applied to equation (A40) or (A4l) for cp or u, and only the contribution 
of the single integral is retained, it is clear that the desired improve- 
ment will not be obtaine d for all Mach numbers up to the critical because 
the function s/l-t^-ku that appears in the denominator of each term 
vanishes with the occurrence of sonic velocity, and the numerator does 
not vanish simultaneously. It is interesting to note, nevertheless, that 
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the result for the pressure coefficient on the airfoil surface that is 
obtained in this way from equation (A4l) , that is. 


Cp(x,0) 




o MJco 

1-lfc, + Cp 


Cp i (x,0 ) 

s/l-M 2 


(A47) 


corresponds to the use of the local, 
rather than free- stream, Mach number 
in the Prandtl-Glauert rule; and that 
this result is the counterpart, in 
transonic small disturbance theory, 
of an approximation proposed by 
Laitone, Szebehely, and Truitt 
(refs. 47 through 51) . It is imme- 
diately apparent that although this 
result differs from the Prandtl- 
Glauert rule in the same direction 
as the higher order approximations, 
the effects of compressibility are 
greatly overestimated at points where 
the local velocity approaches sonic 

velocity. A typical set of results w ^ 

illustrating this statement is shown Moo 

in sketch (n* ) in which are repeated Sketch (n* ) 

the curves of sketch (h) showing the 
variation of Cp with M*. at the midpoint of a 10-percent-thick circular- 
arc airfoil together with the corresponding curve calculated by use of 
equation (A47) • 

The relations developed in the discussion of the elliptic case in 
the main text do not encounter any such difficulties as the local Mach 
number approaches unity. The difference in behavior is a consequence of 
the fact that the latter results are based on the equation for du/dx 
rather than that for cp or u. That such a difference might occur can be 
seen by examin ation of equation (A42) . The denominator of each term again 
approaches zero as the local Mach number approaches 1, but the numerator 
of each term is always zero at the point of maximum velocity. Since sonic 
velocity is first encountered with increasing Mach number at the point of 
mflvirrmm velocity, it is apparent that an indeterminate form occurs at the 
critical Mach number and the possibility at least exists that the gain 
sought by forcing fg to be zero at the point where 5 cr-g/cbc 2 is infi- 
nite will be realized. That a gain in accuracy, although not the mathe- 
matical order, of the solution is actually attained by this procedure is 
shown in the main text by comparison with existing higher approximations. 
Further confirmation of this conclusion is shown by the comparison • 
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illustrated in sketch (o 1 ) in which the 
curves of sketch (h) showing the vari- 
ation with 4 of Cp at the midpoint 
of a 10-percent-thick circular-arc air- 
foil are repeated together with the 
curves calculated "by use of equa- 
tions (a 45) and (A46) representing the 
first two approximations furnished by 
the quadratic method of successive 
approximations . 


Parabolic Case 


The third problem to be discussed is the approximate solution of the 
integral equations for the parabolic case under the restriction that the 
free- stream Mach number is near unity. Two integral equations are given 
depending on whether Ap is positive or negative. Inasmuch as no use is 
made in the main text of the equation for negative Ap, the following 
remarks will be confined to the case for positive Ap/ for which the 
integral equation is given by equation (A20) : 


<P(x ’ Z> - 'wlo ^ § 41 ' ^ JJ ° pfpa5p ' 
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The thi rd, term of equation (Ak8) represents the contribution of 
sources distributed along the surface of any part of the shock system 
that is situated upstream of the point x,z. This term, has the property 
of effectively continuing the source distribution of the first term 
smoothly through a concave comer when the adjacent flow is supersonic. 

In this way, the singularities in the velocity and pressure that occur 
at such comers when the adjacent flow is subsonic are replaced, when the 
adjacent flow is supersonic, by the discontinuous, but finite, jump asso- 
ciated with an oblique shock wave. The contribution of this term vanishes 
if no part of the shock system is situated upstream of x,z or if the 
shock wave is par all el to the z axis. In that which follows attention 
is confined to cases in which it is presumed that one or both of these 
conditions are satisfied for all points situated upstream of the trailing 
edge. The integral over the shock waves thus contributes nothing to cp 
at any point on the airfoil surface, and the remainder of the discussion 
will proceed with considerations involving only the first two terms of 
equation (A48) , and with the corresponding equation for u obtained 
therefrom by differentiation with respect to x. These two equations 
reduce to the following forms for points on the airfoil surface 
(i.e., z=0) : 


cp( x , 0 ) = - — §£= £ X d£ - ~ jj OpfpdRp (M9) 
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It can be seen from examination of the preceding relations that the 
integral equations for the parabolic case are nonlinear and singular just 
as are the integral equations for the elliptic case. The following dis- 
cussion of the approximate solution of the equations for the parabolic 
case will proceed, therefore, through applications of considerations that 
are very similar to those described in the discussion of the elliptic 
case in the preceding section. 

The results found by application of the linearized theory for sonic 
flow described in references 23 through 28 follow from equation (A^9) 
or (A50) by equating fp to zero so that 


u(x,0) = - 


Uco r x az/d| 

^/rtAp uo n/x-£ 


(A51) 


nnri selecting a value for Ap. Various means have been proposed for the 
selection of an appropriate value for. Ap. In reference 28, the only one 
of the above references that pertains directly to two-dimensional flows, 
Maeder and Thommen suggest that Ap be determined by equating it to the 
value for kSu/Sx, obtained by differentiating equation (A51) , that 
occurs at the point along the chord at which u is a m a ximu m in incom- 
pressible flow. As noted in the main text in the discussion of the solu- 
tion for the wedge, the results obtained by application of this linearized 
theory for sonic flow past thin airfoils may be at considerable variance 
with other theoretical and experimental results. A further illustration 
of this statement is provided in sketch (p 1 ) in which the results given 



in sketches (p) , (r) , and (t) for the circular-arc airfoil and the two 
related airfoils that have the point of maximum thickness at 30- and 
70-percent chord are repeated together with the corresponding results 
obtained by application of the procedures described in reference 28. It 
can be seen upon comparison of these results with the experimental results 
shown in sketches (p), (r) , and (t) that the agreement . between the pres- 
sure distribution calculated by application of the linearized theory for 
sonic flow and that measured experl mentally deteriorates as the pressure 
gradient departs from a constant. It is apparent from equation (A51) 
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that the accuracy cannot he improved in any essential ma nn er hy the adop- 
tion of a different rule for the selection of an appropriate value 
for Ap. This follows from, the fact that the entire curve representing 
the pressure distribution is proportional to l/ -J Ap and can he altered 
in scale, hut not in form, hy use of other values for Ap. Possibilities 
for iteration always exist, and it is conceivable that improvements in 
accuracy could he attained hy inserting the solution of linearized sonic 
flow theory into the terms involving fp in equation (Ah9) or (A 50 ) to 
obtain a second approximation, etc. To do so would he a laborious task 
and there is always present a grave danger that the process will diverge, 
or not converge sufficiently rapidly to he useful, when the first approx- 
imation is as far from the proper solution as may he in f erred from 
sketch (p* ) for the airfoil with maximum, thickness at 30 -percent chord. 

The procedures employed in the main text for the approximate solution 
of the equations for the parabolic case closely parallel those used for 
the approximate solution of the equations for the hyperbolic and elliptic 
cases. It is, consequently, not surprizing that the following discussion 
of these procedures from the point of view of the integral equations given 
in equations (A49) and (A 50 ) is very similar to that in the preceding 
sections of the Appen dix . The general considerations are the same in all 
three cases, hut the parabolic case more closely resembles the elliptic 
case than the hyperbolic case because of the singular nature of the kernel 
in the double integral] that is, Op and approach infinity at the 
point |=x, £=z. The expressions applied in the mai n text; follow from 
consideration of equation (A 50 ) as a typical member of an infinite set of 
relations that result if Ap is replaced with the local value of kSu/Bx 
and the contribution of the double integral is disregarded. At free- 
strearn Mach number 1, the function fp thus reduces to zero, at the point 
where Op is infinite, and it again appears plausible that less loss in 
accuracy is incurred by disregarding the contribution of the double inte- 
gral than in alternative procedures in which fp is not zero at this 
point. At free- stream Mach numbers different from unity, fp is not zero 
although it can be made as small as desired by approaching sufficiently 
closely to free-stream Mach number 1. The results obtained by solution 
of the remaining relation, which is a first-order nonlinear ordinary 
differential equation, are completely consistent with the above remarks. 
The pressure distributions calculated by use of only the first approxi- 
mation are indeed in good accord with existing theoretical and experi- 
mental results, and the initial variation of Cp with M^ at free-stream 
Mach number 1 is given exactly; but no indication is provided of the sub- 
sequent variation of Cp with ^ at Mach numbers considerably removed 
from unity. 

If, on the other hand, the pressure distribution is calculated by 
substitution of kcp^ for Ap in equation (A49) and differentiation of 
the resulting expression to obtain an equation for u, the function in 
the double integral that corresponds to fp in the above discussion does 
not vanish at the point |=x, £= z where the kernel is infinite. If the 
preceding discussion can be' considered to imply that the effective 
removal of the singularity is important in the approximate solution of 
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singular integral equations, it may be anticipated that the results 
obtained using equation (A^9) will not, in general, be so good as those 
obtained following the procedure employed in the main text. 
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